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ABSTRACT 

A  set  of  initial  conditions  to  serve  as  primary  inputs  to  the  Department  of 
Defense  Land  Fallout  Prediction  System  (DELFIC)  <s  derived.  The  set  con¬ 
sists  of  lower-boundary  conditions  for  use  by  a  cloud  rise  and  growth  simula¬ 
tion  model.  The  conditions  are  time,  temperature,  soil  burden,  fraction  of  the 
soil  burden  in  the  vapor  phase,  and  size-frequency  distribution  of  fallout  particles. 
The  DOD  Land  Fallout  Prediction  System  predicts  local  fallout  patterns  from  land- 
surface  nuclear  detonations. 
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INTRODUCTION 


The  DELFIC  system  calculations  begin  with  a  simulation  of  the  nuclear  cloud's 
rise  and  growth.  In  this  simulation,  a  set  of  differential  equation,,  which  collect¬ 
ively  describe  the  range  of  change  of  the  system,  Is  solved  over  the  time  interval  of 
loud  rise  history.  These  differential  equations  are  solved  as  an  initial  value 
problem,  1.  e. ,  the  integrations  are  done  numerically  for  a  continuous  succession 
of  small  time  steps  that  cover  the  time  range  of  the  cloud  rise.  Obviously,  a  set 
of  initial  conditions  flower-boundary  conditions)  Is  required  to  begin  the  integrations. 
The  particular  set  of  initial  conditions  discussed  here  is  that  required  by  the  Huebsch 
cloud-rise  model  which  is  described  In  Volume  m  of  this  documentation  and  also  in 
Refa.  1  and  2.  These  conditions  are: 

1.  The  time,  relative  to  detonation  time,  of  the  initial-conditions 
specifications 

2.  The  average  temperature  of  the  nuclear  cloud  at  the  time  of 
the  Initial -conditions  specifications 

3.  The  total  mass  of  soil  material  entrained  in  the  cloud  and  the 
fractions  in  vapor  and  condensed  phases 

4.  Size-frequency  distribution  of  the  condensed  phase  soil  material 

The  initial-conditions  specification  can  account  for  detonation  yields  ranging  from 
0. 01  KT  to  100  MT,  and  heights  of  burst  ranging  from  1  fireball  radius  above 
ground  to  20  scaled  feet  below  ground.  Two  soil  types  are  considered:  siliceous 
and  calcareous.  Implicit  in  the  specification  of  the  Initial  conditions  is  the  realiza¬ 
tion  that  entrained  soil  material  provides  essentially  all  of  the  mass  of  fallout  de¬ 
posited  locally.  This  material  is  made  radioactive  principally  by  occlusion  of 
weapon  debris  material. 

Some  additional  details  of  the  determination,  of  these  Initial  condition,  are 
preaented  in  Ref.  3  which  is  the  final  report  of  a  program  carried  out  at  Tech/Ops 
to  provide  the  initial-condition,  data  set.  The  result,  of  other  work  done  also  at 
Tech/Op.  -  specifically  in  revising  the  initial  temperatures  and  In  determining 
the  fraction  of  the  soil  burden  In  the  vapor  state  -  are  repotted  here  only.  We 
should  mention  that  an  additional  set  of  Initial  conditions  obtained  from  observed 

cloud-rise  data,  such  as  Initial  cloud  dimension,  and  velocity,  Is  to  be  found  In 
Refs.  4  and  5. 
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DETERMINATION  OF  THE  INITIAL  CONDITIONS 

To  supply  the  initial  conditions  by  means  of  a  theoretical  analysis  starting  from 
first  principles  would  require  an  effort  far  beyond  both  the  fiscal  and  temporal  scope 
of  the  program  funded  for  this  work.  In  addition  to  the  need  for  a  detailed  theoretical 
description  of  the  complex  events  and  interactions  during  and  following  a  nuclear  ex¬ 
plosion,  a  very  large  scale  computing  effort  would  be  required  to  provide  a  suffic¬ 
iently  complete  time  history  of  the  fireball  and  early  cloud  development.  The 
analysis  would  need  to  cope  with  a  two-phase  environment  (air  and  ground)  in  such 
fine  detail  and  over  such  a  long  period  of  time  as  to  transcend  the  current  state-of- 
the-art  in  the  field.  Because  of  the  physical  complexity  and  large  scale  of  a  nuclear 
detonation,  it  is  doubtful  that  less  exacting  theoretical  calculations  could  produce 
meaningful  insults.  Attempts  have  been  made  in  the  past  to  provide  the  required 
information  using  simplified  theoretical  models,  but  in  the  end  so  many  questions 
were  left  unanswered  that  one  must  be  extremely  skeptical  about  the  accuracy  of 
the  results.  We  concluded  that  our  best  recourse  was  to  make  use  of  field  data 
observed  at  nuclear  t^sts. 

Needless  to  say  utilization  of  field  data  does  not  solve  all  our  problems.  In 
fact,  the  problems  are  so  numerous  and  involved  that  a  thorough  discussion  of 
them  would  not  be  possible  within  the  space  requirements  of  this  report.  Many  of 
these  are  discussed  briefly  in  the  sections  that  follow;  a  few  of  the  larger  problems 
that  are  only  indirectly  treated  later  are  as  follows.  There  has  not  been  a  really 
large  number  of  nuclear  detonations,  and  detonations  always  have  occurred  under 
restrictive  conditions.  For  example,  U.  S.  test  shots  of  interest  to  fallout  studies  have 
occurred  in  only  two  environments:  the  Nevada  Test  Site  and  the  Pacific  Testing 
Grounds.  The  former  has  a  desert  environment;  the  latter,  a  tropical  ocean  en¬ 
vironment.  All  of  the  large-yield  shots  have  taken  place  in  the  Pacific  whereas 
most  of  the  small-yield  shots  have  occurred  in  Nevada.  Most  of  the  Nevada  shots 
have  been  air  bursts  or  tower  bursts  which  produced  little  or  no  local  fallout.  At 
only  a  few  shots  have  fallout-study  programs  been  carried  out  on  a  scale  large 
enough  to  yield  comprehensive  results;  several  of  these  studies  were  done  in  the 
Pacific  in  connection  with  multi-megaton  detonations  that  scattered  fallout  over 
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•-uch  large  areas  of  ocean  that  only  sparae  samplings  could  be  taken.  There  are 
essentially  no  applicable  data  for  the  lowest  end  of  the  yield  range,  and,  of  course 
no  known  shots  have  been  as  powerful  as  those  at  the  highest  end. 

On  the  other  side  of  the  coin,  we  have  been  fortunate  in  having  available 
several  sets  of  extremely  useful  data.  Particularly  noteworthy  are  the  unpublished 
cloud  rise  and  growth  data  of  Edgerton,  Germeshausen  and  Grier,  Inc  6  the 
temperature  data  of  HUlendahl,  the  crater-scaling  equations  developed  byNordyke 

a  allout  data  collected  from  shots  TEAPOT  Ess9  and  SUN  BEAM  Small  Boy,  10and 
the  DASA  compilation  of  fallout  data. 11 

Time  of  Initial  Condi  dons  Specifications 

The  first  step  in  establishing  a  set  of  initial  conditions  is  to  select  an  appropriate 
me  at  whic  i  the  initial  conditions  are  tc  be  specified.  Varic-s  alternatives  are 
available.  For  example,  we  could  select  the  detonation  time  but  then  we  would 
necessarily  have  to  include  In  our  subsequent  analyses  and  computational  simula¬ 
tion.  the  very  early  phenomenology  of  the  radiative  fireball  growth  in  audition  to 
cloud  rise  and  growth.  A  second  alternative  would  be  to  select  the  time  at  which  a 
particular  temperature  Is  reached.  TOs  alternative  is  particularly  attractive  be¬ 
cause  one  might  pick  a  temperature  at  which  ft  is  certain  that  all  of  the  entrained 
.0,1  material  has  condensed.  Obviously,  this  simplifies  the  ensuing  cloud-dynamics 
calculations  and,  even  better,  it  eliminate,  the  troublesome  problem  of  establishing 
e  faction  of  the  soil  burden  in  the  vapor  state.  Another  convenient  temperature 
might  be  one  that  is  critical  in  the  partitioning  of  activity  between  refractor,  and 
volatile  radioactive  Isotopes.  This  would  link  the  initial-time  determination  with 
the  activity  fractionation  process.  Actually,  we  have  used  none  of  these  alternatives. 

to  choosing  a  time  for  the  Initial-conditions  specification,  one  should  not  lose 
sight  of  the  primary  purpose  lu  supplying  the  initial  conditions:  namely,  ,o  provid- 
a  set  of  lower-boundary  condition,  for  use  by  a  dynamic  cloud  rise  and  growth 
simulation.  With  this  purpose  In  mind,  a  study  of  available  cloud  rise  and  growth 
data  derived  from  movie  films  of  test  shot.6  readily  reveals  an  appropriate  time 
However,  to  thoroughly  understand  the  significance  of  the  time  selected,  a  brief 
discussion  of  the  cloud  data  ,n  terms  of  the  dynamic,  of  cloud  rise  and  growth  Is 
warranted.  Figure  1  Illustrate,  typical  plots  of  cloud  dimensions  versus  time. 
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We  have  found  it  convenient  and  appropriate  to  divide  the  development  of  a  nuclear 
cloud  into  four  phaaee.  These  are  (in  temporal  sequence): 

1.  A  fireball  phase 

2.  A  pseudohydrostatic  cloud  rise  phase 

3.  A  similarity  cloud  rise  phase 

4.  A  final  expansion  phase  leading  to  stabilization. 

The  fireball  phase  covers  the  interval  from  the  detonation  time  until  radiative 
growth  of  the  fireball  has  virtually  ceased  and  pressure  equilibrium  with  the  atmo¬ 
sphere  surrounding  the  fireball  has  been  achieved.  The  prompt  effects,  including 
the  electromagnetic  pulses  and  the  close-in  blast,  have  largely  subsided  by  the  end 
of  this  period.  During  phase  2,  the  fireball  is  still  hot  (in  fact,  the  second  temper¬ 
ature  maximum  usually  occurs  early  in  this  phase),  but  not  hot  enough  for  radiation 
to  dominate  the  cooling  process.  Graphs  of  fireball  center  height  versus  time  indi¬ 
cate  decelerative  rise.  Because  of  rigorous  vortex  circulation,  the  fireball  entrains 
appreciable  amounts  of  ambient  air;  this  cause*  rapid  cooling.  The  entrainment  rate 
during  phase  2,  though  significant,  is  substantially  less  than  the  rate  that  occurs 
subsequently  and,  in  comparison  with  this  later  behavior  (phase  3),  the  cloud  be- 

more  M  bubble  rising  in  water.  ^  For  this  reason  we  have  labeled  it 
the  pseudohydrostatic  rise  phase.  The  comparatively  low  entrainment  rate 
probably  is  caused  by  the  fact  that  the  cloud  temperature  is  sufficiently  high,  and 
as  a  consequence,  the  difference  in  density  between  the  cloud  gases  and  the  ambient 
Ur  1.  nfflctntly  *r«t  (punhlM„/pcloud  >  10),  to  prevent  an  effective  momentum 
transfer  from  die  circulating  cloud  gases  to  the  ambient  air.  An  interpretation  of 
this  type  has  been  presented  previously  by  Batchelor. 13 

A  short  period  ot  accelerative  cloud  center  rise  appears  to  occur  in  a  small 
time  interval  centered  on  the  boundary  between  phases  2  and  3.  Subsequent  to  this 
period,  during  phase  3,  decelerative  rise  and  growth  is  resumed.  An  appreciation 
of  the  critical  influence  of  vortex  circulation  on  cloud  rise  and  growth  is  required 
for  an  understanding  of  these  processes.  Figure  2  shows  an  artist’s  conception  of 
the  geometry  and  structure  of  organised  vortex  circulation  within  a  nuclear  cloud. 
(Detailed  discussions  of  the  physics  of  vortex  rings  and  the  efficacy  with  which  vortex 
ring  geometry  and  circulation  fit  observed  nuclear  cloud  behavior  have  been  presented 
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Figure  2  Toroidal  Circulation  Within  a  Nuclear  Cloud  (From  Ref.  14) 

by  Norment. 15 * 16  Movie  films  of  nuclear  clouds  show  that  during  phase  2  this 
circulation  is  very  vigorous  but,  as  noted  earlier,  because  of  the  gross  density 
difference  between  the  cloud  and  ambient  gases,  there  is  little  effective  transfer 
<rf  circulatory  momentum  from  the  cloud  to  the  surrounding  air.  However,  at  approxi 
mately  the  time  of  the  beginning  of  phase  3  the  combined  influence  of  the  lowered 
temperature  and  the  circulation  vigor  are  sufficient  to  cause  more  effective  momen¬ 
tum  transfer.  When  this  occurs,  the  vortex  circulation  causes  increasingly  vigorous 
sxpunsion  of  the  cloud  into  the  surrounding  air.  This  dynamic  interaction,  which 
causes  much  higher  rates  of  cloud  growth,  continues  throughout  phase  3.  The  rise 
«ad  growth  are  correlated  during  phase  3  in  such  a  manner  that  geometric  similarity 
of  cloud  shape  is  maintained  approximately.  That  is  to  say,  if  we  view  the  cloud  as 
an  oblate  ellipsoid,  the  ellipsoid  eccentricity  does  not  vary  greatly  during  phase  3. 

Subsequent  to  phase  3,  cloud  behavior  1.  dependent  on  weapon  yield  and  height 
erf  burst.  For  small  detonations  at  low  altitudes,  rise  momentum  and  circulation 
•ubeide  at  a  level  below  the  tropopause.  For  Urge  detonations,  the  tropopause 
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where  W  is  In  Idlotons.  No  height  of  burst  effect  is 
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apparent  in  the  data. 
°an  be  computed  from  the  equations  of  Hlllendahl. 8 
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Hons  con  be  combined  by  linearly  Interpolating  between  them  to  yield 
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where  X  is  the  scaled  height  of  burst  in  units  of  ft/KT1^3. 

^  °f  ^  <4)  ,0r  sub,u  -»  bural»  >»«  not  been  adeqna  ely  Justified  since 
appropriate  cloud-ria.  an  j  ^  date  are  no.  available  for  subsurface  bursts. 

Figure  3  gives  computed  values  of  t(  as  a  function  of  yield  and  ;  burst 
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Figure  3.  Initial  Tim.  vs.  Yield  a.  a  Function  of  Scald  Height  of  Burst 
The  Initial  Temperature 

After  a  particular  Urn.  a.  which  to  specify  the  initial  conditions  ha.  been 
ected,  our  next  task  is  to  estimate  an  svenm  * 

.  „  “ ln  avei*g©  cloud  temperature  at  that  time  We 

have  available  both  theoretical  means  and  observed  data  hv  wM  k  ♦ 

^  ooserved  data  by  which  to  accomplish  this. 

The  problem  cm,  be  attacked  theoretically  at  either  of  two  widely  separated  level,  of 

^  A‘ <™-«-....mp.ep.Lon,r-C. 
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energy  yield  between  blast  and  shock,  prompt  radiation,  nuclear  energy,  and  energy 
used  in  heating  and  lifting  the  fireball.  The  cloud  temperature  at  the  time  of  inter¬ 
est  is  then  determined  from  the  residual  energy  after  accounting  for  all  the  losses. 

17 

Moulton  has  pointed  out,  however,  that  because  of  complex  interactions  between 
these  processes,  simple  energy  partitioning  is  not  physically  meaningful.  Further¬ 
more,  even  if  the  partitioning  were  possible  in  principle,  inevitably  there  would  be 
large  uncertainties  in  each  term  of  the  energy-balance  equation.  Since  the  information 
we  seek  is  obtained  from  differences  between  energy  terms  witn  largo  uncertainties, 
we  can  expect,  in  many  cases  at  least,  that  our  results  will  be  of  the  same  order  as 
the  uncertainties.  Such  results,  then,  could  not  be  considered  reliable.  The 
theoretical  approach  at  the  higher  level  of  sophistication  requires  numerical  solu¬ 
tions  of  a  set  of  basic  hydrodynauiic  and  energy-balance  differential  equations  that 
can  describe  the  processes  in  action  from  detonation  time  to  the  time  of  the  initial  - 
conditions  specification. .  The  model  used  would  need  to  cope  with  the  air-ground 

interface  and  consequently  would  have  to  treat  a  two-material  environment.  This 

18 

type  of  calculation  is  grossly  demanding  of  the  most  advanced  computers  and 
would  be  well  beyond  both  the  temporal  and  fiscal  scope  of  this  effort.  Fortunately, 
there  are  some  excellent  experimental  temperature  data  observed  from  fireballs  ol 
nuclear  test  shots,  and  our  best  course  is  to  utilize  these  data. 


The  experimental  data  used  are  Hillendahl's  power  temperatures  listed  in 
Ref.  7.  These  power  temperatures  were  derived  by  Hillecdahl  from  observed 
calorimetric  data  by  using  the  Stefan-Boltzmann  power  emission  law  in  the  form 


where 


T  = 


_p2 

*R  ST 


1/4 


(5) 


T  *  the  power  temperature 

12  -2  -1  4 

<t  -  Stefan's  constant  (1.356  x  10  cal  cm  *ec  deg) 

R  *  the  slant  range 

-2 

Aq  »  the  measured  energy  in  the  time  interval  of  length  At  (cal  cm  ) 
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0  =  the  wel«hted  cross-sectional  area  of  the  fireball  observed  (cm2) 

*  ■  a  transmission  factor  (0  <  $<  l) 

€  *  average  emissivity  of  the  fireball  surface  (taken  as  1). 

The  power  temperatures  are  given  as  functions  of  scaled  time,  J. e. ,  t/t  ,  where 
l2m  18  1116  tlme  of  the  second  temperature  maximum.  2m 

Plots  of  log  T  versus  logt/t^  almost  always  are  linear  in  the  range  of 
scaled  time  between  about  1. 1  through  4. 0.  At  larger  values  of  scaled  time,  the 

magnitude  of  the  slope  often  begins  to  decrease.  We  have  found  that  individual  sets 
of  data  can  be  fitted  to  an  equation  of  the  form 


T  "  K  +  1500  . 

Xl2m/ 


A  crude  theoretical  justification  for  Eq.  (6)  is  given  in  Apoendix  A.  An  illustrative 
graph  of  observed  and  calculated  temperatures  is  shown  in  Figure  4. 
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Up  to  this  point,  the  work  reported  here  is  the  same  as  that  given  in  Ref.  3; 
however,  since  its  publication,  additional  work  has  been  done  that  has  led  to  a  re¬ 
vision  of  the  temperature  equations  and  curves  reported  in  Ref.  3.  This  additional 
work  was  made  possible  by  an  independent  program  of  studies  at  Tech/Ops  of  ob¬ 
served  kinematic  cloud-rk  j  data  that  resulted  in  determinations  of  cloud  volumes 
at  our  initial  time  as  a  function  of  yield  and  height  of  burst.  *  With  this  new  informa¬ 
tion  available  we  were  able  to  perform  energy  balance  calculations  which  indicated 
that  our  temperatures  for  low  yield  surface  and  subsurface  burst  were  too  high.  In 

light  of  this  information,  we  have  reconsidered  our  analysis  of  Hillendahl's  tempera¬ 
ture  da  o’,  ns  follows. 

In  our  original  study,  the  constants  K  and  n  were  determined  graphically  for 
each  of  a  number  of  shots.  To  determine  yield  scaling  for  theae  parameters  that  is 
independent  of  height  of  burst  effects,  it  was  necessary  to  separate  the  data  into  two 
groups:  one  for  air  bursts  and  the  other  for  surface  bursts.  Only  for  these  two 
categorizations  are  the  data  sufficiently  numerous  to  permit  scaling.  Therefore  n 
and  K  were  determined  as  function,  of  yield  separately  for  the  air  bursts  and  sur¬ 
face  bursts.  To  determine  height  of  burst  between  these  extremes,  a  linear  height 

of  burst  interpolation  was  assumed  arbitrarily  for  both  K  and  n.  (Additional  details 
are  given  in  Appendix  A  of  Ref.  3. ) 

In  determining  yield  scaling  of  n  and  K,  we  observed  that  the  water  surface 
burst  data  alone  yielded  relatively  clean  linear  log-log  graphs  whereas  data  from 
land  surface  (which  in  all  cases  included  considerable  water  also)  and  shielded  bursts, 
while  not  substantially  different  from  water  surface  data,  tended  mostly  to  increase 
the  scatter  on  the  graphs  and  thereby  add  confusion  to  a  situation  that  already  was 
somewhat  obscure.  Because  of  this  we  used  only  the  water  surface  data.  In  rework¬ 
ing  the  data,  we  included  these  "wet"  land  surface  and  shielded  detonations.  These 
amounted  to  six  additional  shots  and  included  one  at  a  much  lower  yield  ^ibout  40  KT) 
than  was  available  in  the  initial  data  set.  The  augmented  data  set  does  uot  display  a 
clearly  defined  height  of  burst  effect  in  the  n  versus  W  graphs  so  that  no  such  effect 

for  n  appears  in  the  final  equation.  The  height  of  burst  effect  for  K  remains  clearly 
evident. 
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The  final  parameters  decided  upon  did  not  account  for  different  types  of  sur¬ 
face  material  since  no  unambiguous  data  are  available.  (All  land-surface  bursts  in 
the  Pacific  incorporated  some  water. )  The  linear  relationship  in  X.  assumed  for  the 
effect  of  height  or  depth  of  burst,  could  not  be  checked  because  insufficient  data 
exist  for  tests  with  scaled  burst  heights  between  surface  and  air  bursts,  and  no 
information  on  temperatures  is  available  for  underground  bursts.  The  range  of 
yields  did  not  extend  to  the  extremes  at  the  high  and  low  ends  of  the  scale  (0. 01  KT 

and  100  MT),  so  these  values  are  the  results  of  simple  linear  extrapolations  from 
available  data. 

The  complete  expression  for  calculating  the  initial  temperature  given  by 
Eq.  (6)  is 

T  -  K  (')  +  1500 
Vl2m  / 

in  which 


and 


where 

K  «  5980. 
o 

k  --  1.145 
pQ  =  -0.0395 
p  -  0. 0264 
nQ  ■  -0.447 
qQ  *  0. 0436  ; 

W  is  the  detonation  yield  (KT)  and  X  is  the  scaled  height  of  burst  <ft/KT1/3). 
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The  height  of  burst  and  yield  ranges  of  the  modei  are 


and 


-20  <  X  <  i80 


0.  01  <  W  <  100,000 


Temperatures  obtained  from  Eqs.  (6),  (7),  and  (8)  are  to  be  used  as  average  cloud 


temperatures.  Figure  5  shows  curves  of  temperature  at  t{  as  a  function  of  yield 
and  height  of  burst. 


worn 


Figure  5.  Calculated  Cloud  Temperature  at  t.  vs.  Yield  as  a 
Function  of  Height  of  Burst1 
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It  Is  evident  that  the  prime  factors  influencing  the  soil  loading  of  the  fireball 
are  the  yield  and  height  (or  depth)  of  burst  of  the  nue’ear  device.  In  addition,  the 
characteristics  of  the  medium  in  which  the  device  is  detonated,  the  topography  of 
tl.e  blast  site,  and  the  nature  and  quantity  of  other  materials  in  the  immediate  en¬ 
vironment  of  the  blast  are  factors  which  also  contribute  substantially  to  the  final 
fallout  patterns. 

There  are  a  number  of  procedures  that  conceivably  can  be  used  to  estimate 
soil  contents  of  nuclear  fireballs.  We  have  chosen  to  label  these  methods  (1)  en- 
engy  partition,  (2)  hydrodynamics  calculation,  (3)  cloud  buoyancy,  (4)  back  extra¬ 
polation,  (6)  specific  activity,  (6)  crater  volume,  and  (7)  fireball  volume.  We  have 
studied  and  evaluated  each  of  the  methods;  discussions  of  the  results  of  our  deli¬ 
berations  are  presented  in  Ref.  3  and,  therefore,  will  not  be  repeated  here.  In 
determining  soil  loadings  for  specific  test  shots  and  also  in  deriving  yield  and 
height  of  burst  scaling  functions  we  have  used  methods  4,  5,  6,  and  7. 

Soil  Burden  Determination «  fox  Selected  Test  Shots.  In  constructing  and 
calibrating  our  model,  we  were  restricted  to  a  consideration  of  nine  shots  whose 
yields  ranged  from  1. 2  to  15, 000  KT,  with  scaled  heights  (depths)  of  burst*  extend¬ 
ing  from  -63.  51  ft  .  (KT)'l/3, 4  to  130. 152  ft  •  (KT)"l/3* 4.  (This  restriction  to 
nine  shots  was  dictated  by  the  extent  of  information  available  regarding  radiological 
and  fallout  mass  collection  data. )  For  two  of  the  nine  shots,  the  available  data  are 
of  Interest  but  are  not  quite  appropriate  for  our  needs.  Specific  activity  data  for  one 
shot  were  taken  from  fallout  samples  collected  on  a  Japanese  fishing  boat;  thrir  ac¬ 
curacy  i*  open  to  question.  Another  case  considered  is  a  barge  shot  over  shallow 
water  which  consequently  cannot  properly  be  classified  as  a  land-surface  burst. 

We  used  two  methods  to  reduce  the  raw  data.  The  first  method  directly  inte¬ 
grates  the  mass  deposited  over  the  local  fallout  area.  Continuous  curves  of  fallout 
mass  per  unit  area,  m,  desalted  at  distances  R  along  radial  lines  emanating  from 
ground  zero  were  obtained  from  fallout  collection-station  data  by  leaot  squares  fit¬ 
ting  of  the  data  to  a  function  of  the  form 

In  the  soil  faurdm  work  we  have  used  the  scaled  height  of  burst,  A,  defined 
by  Nordyke:8  A  »  hW-l/8. 4,  where  h  is  the  height  of  burst  in  feet  and  W  is  yield  in 
kilotons.  Note  that  for  the  time  and  temperature  determinations  we  have  used  the 
scaled  height  X=  hW1/3. 
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!n(m)  =  afn(R)+b  , 


(9) 


where  a  and  b  are  conatanta  that  are  unique  to  a  particular  radial  line.  Using  this 
function  for  extrapolation,  we  were  able  to  construct  lsomass  contours  well  beyond 
the  actual  range  of  observation,  and  thereby  extend  the  mass  integrations  to  include 
virtually  all  of  the  local  fallout.  The  derivation  of  Eq.  (9)  and  a  description  of  the 
data- reduction  method  are  given  in  Appendix  B  of  Ref.  3.  This  method  was  used  for 
two  of  the  nine  shots  rince  only  for  these  two  were  a  sufficient  number  of  collection 
stations  employed.  These  shots  are  SUN  BEAM  Small  Boy  and  TEAPOT  Ess. 

Figure  8  gives  the  extended  isomass  contours  for  TEAPOT  Ess  (1.3  KT,  -67  ft). 


FlgUr*  C01*"1"  CThe  contour  levels  are 

1000,  100,  50,  10,  1.0  ,  0.5,  0.1  g/ft2 
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For  three  REDWING  shot*  -  Tew*,  Zuni,  and  Lacroese  -  total  soil  masses 
in  the  clouds  were  obtained  by  analysis  of  specific  activity  data*  that  were  measured 
from  fallout  samples  collected  for  these  shots.  Specific  activities  for  Ce144  were 
used  instead  of  Mo  data  for  reasons  discussed  in  Section  4. 9  of  Ref.  20.  Selected 
data  recommended  by  Mr.  R.  C  Tompkins,  Nuclear  Defense  Laboratory,  from 
Tables  3. 21,  3. 24,  and  3. 25  erf  Ref,  20  were  used.  Specific  activities,  in  units  of 
flssions/unit  mass  of  sample,  for  the  remaining  four  shots  were  taken  from  Miller21 
(TaUr  3. 11);  these  specific  activities  were  based  on  Mo"  data. 

For  SUN  BEAM  Small  Boy,  radiological  data  as  well  as  mass  deposition  data 
were  obtained  at  the  fallout  collection  stations.  This  afforded  the  opportunity  to 
compare  masses  estimated  by  the  two  methods.  If  we  exlude  the  results  from  two 
stations  which  appear  to  be  spurious,  the  average  mass  determined  by  the  specific 
activity  method  agrees  with  the  value  obtained  by  integration  of  the  'all out  mass 
deposition  data  to  within  40%. 


Major  sources  of  error  in  the  two  methods  of  data  reduction  are  known,  but 
quantitative  determinations  of  the  errors  are  not  possible  since  the  total  fallout 
from  a  nuclear  explosion  has  never  been  determined.  In  the  method  of  direct  mr.ss 
integration,  the  largest  error  lies  in  the  assumption  that  the  mass  measured  on  the 
*  ound  contains  all  the  fallout,  and  thus  all  the  soil  mass  in  the  cloud.  A  sizable 
proportion  of  the  lightest  particles  in  the  cloud  settle  well  beyond  the  range  of  the 
cloee-in  fallout  area.  In  fact,  many  particles  are  carried  far  aloft  and  contribute 
to  worldwide  fallout.  This  effect  is  particularly  enhanced  if  the  winds  above  the 
blast  site  are  strong.  Errors  of  this  type  are  expected  to  be  minimized  for  under¬ 
ground  shots  and  shots  in  soil  with  relatively  Urge  soil  particle  sizes,  since  we  then 
expect  lower  cloud  heights  and  faster  particle  fall  rates.  A  second  major  source  of 
error  is  caused  by  contamination  in  the  collectors.  Since  fallout  areas  are  too  hot 
to  be  entered  immediately,  ambient  winds  naturally  will  add  dust  to  some  collectors 
*****  the  samples  can  be  recovered.  For  SUN  BEAM  Small  Boy,  which  was 


l m<Kh0d  COnalttM  <*  extrapoUttng  measured 

^ >  conversion  from  curies/mg  of  sample  to  equivalent 
deteMtton!^  tIrl,l0n  tot0  to*!  number  of  fissions  produced  by  the 


17 


eroded  over  tell  in  the  telc-.lxe  tang.,  both  source,  of  error  may  hove  «riou.iy 
bta«d  the  obeerved  .ire  di.tribution,.  Although  TEAPOT  Ess  *..  detonated  in  M„ 
which  had  larger  granule.  than  SUN  BEAM  Small  Boy,  there  was  a  Use  surge  which 
extended  out  well  over  a  mile  from  ground  aero.  The  c.oae-in  station,  therefore  may 

*  m,terial  not  properly  be  classified 

The  main  .ouree  of  error  in  determining  soil  mas.  loading  of  the  cloud  by 
activity  meamirement.  lie.  in  determining  an  average  value  of  the  specific  activity. 
Specific  activity  will  vary  with  particle  .ire  and  location.  Since  tie  number  of  col- 
ectlon  .tattoo,  for  the  test,  for  which  thl.  method  wa.  used  were  too  few  to  ac¬ 
curately  determine  the  average  .pacific  activity  with  respect  to  location,  we  were 

forced  to  amune  that  the  measurement,  actually  available  were  representative  of 
the  avenge  apecific  activities. 

Pgvelopment  of  Scaling  Function..  After  the  total  fallout  mas.  had  been  de¬ 
termined  for  the  nine  case.,  tern  model.,  one  for  ,„b.urface  and  the  other  for  air 

’  <teV8l°ped  10  correl*te  <*>*  PPmSte  with  yield  and  height  of  burst  We 
necessarily  relied  imavily  on  ob.erved  phenomenology  associated  with  nuclear  blast. 

or  .urfwe  and  .ub.urf«.  bur.te  ..  found  It  natural  to  ask  whether  mass  lifted 
could  be  correlated  with  crater  voiume.  We  plotted  fallout  mas.  verso,  crater 
vohmm.  -calculated  for  a  paraboloid-shaped  crater  (from  crater  data  give.  In 

”*\  ,k  ™**nph  •h0W,  11  '*  <0  suppose  that  a  power-law  rela¬ 

tionship  hold,  between  the  ma..  lifted  mm  the  volume  of  the  resulting  crater  The 

•lope  of  the  learn  wpure.  .tralght  line  fitted  to  the  date  1.  0. 86.  since  the  plot  I. 

«.  a  tog-log  .cm.,  tel.  me,  implies  teat  to  the  extent  that  our  date  are  reliable 
w.  nuty  Mtert  that  the  ma..  lifted  in  a  detonation  I.  approximately  directly  pro-' 
porttaual  to  the  volume  cf  the  crater  produced. 

Tlte  problem  of  the  yield  and  ttopth  of  bum  dependence  of  the  soil  loading  1. 

"T  rr  ,lDC'  ”  °Md  l00k  mly  *l  U“  of  crater  volume. 

than  tat^T.l  l"r“'  “  "r“  **  h“  b”"  "”"*“*»tP‘t  much  more  thoroughly 

8PeC‘,IC*Uy'  U*  -  c cater  volume 

to  yield  and  depth  cf  bum  ha.  been  InveWgatm!  In  detail  by  Nordyke. 8  While  the 

hwctKrn.  he  derive,  are  hated  only  on  high  explclve  bur.te,  hi.  comp. ri.cn  of  the 
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function*  with  data  from  low-yield  nuclear  bursts  shows  compatibility,  the  dis¬ 
crepancies  in  crater  dimensions  being  no  more  than  10%.  One  obvious  difficulty 
in  calibrating  a  volume  scaling  function  is  die  lack  of  really  pertinent  data  at  very 
high  yields.  Because  of  this,  and  since  data  on  craters  produce  j  by  high-yield  ex¬ 
plosions  are  in  general  sparse,  we  have  chosen  to  use  Nordyke's  scaling  functions 
extrapolated  to  the  high-yield  region.  Nordyke's  functions  correlate  quite  well 
with  the  available  crater  data  for  the  high  -yield  region.  We  have,  then: 

M  =  K  R2(z)D(z) W3^3, 4  ,  (10) 

z 

where 

M  =  the  total  mass  of  soil  in  the  cloud  (g) 

W  *  the  yield  (KT) 

z  *  the  scaled  depth  of  burst  0ft  •  KT-1^3, 4) 

R(z)  -  112. 5  +  0.  755z  -  9. 6  x  10~  6  z3  -  9. 1 1  x  l(f  12z5 
D(z)  -  32.  7  +  0.  851z  -  2. 52  x  l(f  5  z3  +  1.  78  x  10_10z5 

and  K  is  a  constant  (discussed  later)  that  relates  the  crater  volumes  to  soil  n»»§s 
z 

lifted  in  the  cloud. 

For  height  of  burst  greater  than  zero,  no  study  has  been  made  of  crater  vol¬ 
umes  comparable  to  that  of  Nordyke's  for  surface  and  subsurface  bursts.  Further¬ 
more,  crater-dimension  data  that  are  available  do  not  correlate  with  yield  in  a 
consistent  manner.  Therefore,  a  different  scaling  model,  based  on  the  following 
rationalizations,  was  used  for  air  bursts.  The  most  Intense  blast  and  thermal 
effects  of  nuclear  detonations  are  contained  within  the  region  of  the  fireball,  and 
it  is  within  the  reach  of  this  region  that  the  cratering  and  soil  incorporation  take 
place.  Of  course,  the  dimensions  of  the  fireball  vary  with  time  and  never  are 
precisely  defined.  However,  we  know  that  the  fireball  has  a  spherical  shape,  and 
we  can  reasonably  assume  that  its  volume  is  approximately  proportional  to  the 
detonation  yield.  Therefore,  a  characteristic  dimension  of  the  fireball,  for 
example  the  radius,  will  vary  approximately  as  the  cube  root  of  the  yield,  namely 
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r  -  W  .  (This  relation  provide*  the  bests  for  the  form  of  the  height  of  burst 
and  depth  of  tMirst-acellng  functions,  >  =  hW'l/3  (Ref.  14.  sections  3. 55-3. 58). 
or  A  hW  (Ref.  8). )  If  the  cratering  and  soil  incorporation  occur  within 
reach  a  the  fireball,  it  is  reasonable  to  expect  that  the  mass  of  soil  lifted  will  be 

proportional  to  the  volume  of  the  fireball  segment  that  appears  to  be  intersected 
by  the  ground  (see  Figure  7). 


—  Ground  Surface 

Figure  7.  Intersection  of  a  Fireball  with  the  Ground 

Our  equation  for  estimating  soil  mass  lifted  by  a  nuclear  cloud  has  the  form 
of  the  equation  for  the  volume  of  a  segment  of  a  sphere.  This  segment  may  be 
thought  of  as  that  portion  of  fireball  volume  that  is  intersected  by  the  ground.  In 
our  equation  the  effects  of  yield  and  height  of  burst  are  expressed  in  the  combined 
form  of  the  scaled  height  of  burst  A.  Figure  7  shows  the  pertinent  geometry  for 
a  fireball  with  radius  r  and  with  its  center  at  an  altitude  h  above  the  ground.  The 

hatched  area  represents  the  region  of  the  fireball  that  intersects  the  ground,  the 
volume  of  which  is  given  by 


V  -  j(r-h)2  (2r  +  h)  . 

To  Opacify  a  particular  nlue  for  r,  and  at  the  aame  time  provide  a  reiationahip 
between  b  and  W  that  Impooea  a  properly  ecaled  limit  on  h  above  which  no  local 

Mlout  can  ooonr,  we  note  (Ref.  14,  Section  2. 118)  that  local  fallout  occura  only 
when  the  relation 


h  <  180W°* 4 


(12) 
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is  approximately  obeyed.  This  relation  provides  us  with  an  appropriate  exponent 

to  use  on  the  yield  in  oar  scaling  function;  however,  to  be  consistent  and  to  ensure 

continuity  at  the  ground  surface  between  the  scaling  functions  for  subsurface  and 

l/3  4 

air  bursts,  we  have  chosen  to  continue  to  use  W  ‘  .  Now,  we  define  a  radius¬ 
scaling  function  by 

r  -  180  W1^3* 4  (13) 

In  terms  of  the  scaled  height  of  burst  and  radius-scaling  functions,  Eq.  (13) 
becomes 

V  *  §  W3^3, 4  (180  -  A)2  (360  +  A)  .  (14) 

Thus  we  obtain  finally 

M  -  Ka  W3/3.4 

q80  -  A)2  (360  +  A)  .  (15) 

The  constants  K  of  Eq.  (10)  and  K.  of  Eq.  (15)  were  determined  as  follows. 

We  noted  that  of  all  detonations  studied,  the  mass  data  for  TEAPOT  Eos  are  by 

far  the  most  reliable.  Therefore,  we  used  this  single  shot  to  calibrate  the  models. 

K  was  determined  by  direct  solution  of  Eq.  (10).  Then,  to  satisfy  the  require- 
z 

ment  of  continuity  at  z  ■  A  ■  0,  we  set  Eqs.  (10)  and  (15)  equal  and  thereby  ob¬ 
tained  K^.  In  light  of  the  relative  unreliability  of  the  remainder  of  the  observed 
data,  this  means  of  calibration  is  at  least  as  valid  as  any  other  that  might  be 
proposed.  Figure  8  shows  calculated  masses  a*  a  function  of  yield  and  height  of 
burst. 

The  main  justification  for  our  models  is  the  agreement  between  the  calcu¬ 
lated  and  the  measured  data.  It  is  Interesting  to  note  that  the  Stanford  Research 
22 

Institute  experimentally  determined  that  masses  of  soil  lifted  were  proportional 
to  crater  volumes  for  a  series  of  low-yield  surface  shots  which  they  examined. 
They  also  determined  that  crater  volume  is  proportional  to  yield  to  the  0. 92  power 
which  is  in  oloee  agreement  with  our  yield  to  the  S/3. 4  «  0. 88  power.  Also,  since 
the  work  reported  here  was  completed,  results  of  a  somewhat  similar  study  at  the 
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JO  -to 


Figure  8.  emulated  Total  Soil  Maas  Lifted  by  the  Cloud  as  a  Function 
of  Yield  for  Various  Scaled  Heights  of  Burst 

23 

Lawrence  Radiation  Laboratory  have  been  released  T  >tal  masses  of  soil  mater¬ 

ial  carried  in  clouds  for  a  number  of  CASTLE  and  REDWING  shots  were  determined 
by  radiochemical  analyses  of  samples  collected  by  aircraft  traverses  through  the 
clouds.  For  the  four  cases  where  these  data  can  be  compared  directly  with  ours, 
mass  ratio  (ours/theirs)  la  0.  7  for  the  best  case  and  27  for  the  worst. 

Hie  complete  set  of  equations  governing  the  soil  burden  model  are  as  follows: 
MSoil  "  kewS/3,  4r2<z>  d<z)’  0  <  z  <  20 

MSotl  *  KaW3/3,4(180-A)2(360+  A),  180>A>0 
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where 


M  *  the  total  mass  of  soil  in  the  cloud  (g) 

W  -  the  yield  (KT) 

A  *  scaled  height  of  burst  (Tt/KT1^3, 4) 
z  -  scaled  depth  of  burst  (ft/KT1^3,4) 

R(z)  «  112. 5  +  (7. 55  x  10"1)  z  -  (9.  6  x  10'6)  z3  -  (9. 11  x  10'12)  z5 
D(z)  -  32.  7  +  (8.  51  x  10*1)  z  -  (2.  52  x  10“6)  z3  +  (1.  78  x  10‘10)  z5 
«  77. 41  g/ft3 
Kz  »  2182  g/ft3. 

Phase  Partitioning  of  the  Soil  Burden.  In  addition  to  the  average  cloud  temp¬ 
erature  and  the  total  soil  burden  of  the  cloud,  certain  other  information  is  essential 
to  properly  account  for  subsequent  cloud  rise  and  growth,  distribution  of  activity 
among  the  particles,  and  the  atmospheric  transport  of  fallout.  The  thermodynamic 
cloud-rise  computations  must  have  the  quantity  of  soil  material  in  the  vapor  phase 
and  the  average  temperature  of  condensed  phase  material.  A  knowledge  of  the  fall¬ 
out  particle  size-frequency  distribution  is  essential  to  both  the  distribution  of 
aotivity  on  the  particles  and  atmospheric  transport  computations.  To  produce  high 
quality  estimates  of  the  quantities  mentioned  would  require  an  extensive  combined 
theoretical  and  computational  effort  that  would  start  from  first  principles  and  fol¬ 
low  the  many  relevant  complex  processes  from  essentially  shot  time  until  the 
early  oloud  development.  Such  an  effort  never  has  been  undertaken. 

These  particular  problems  are  complicated  by  several  considerations.  At 
the  time  of  Interest,  tj,  there  is  every  reason  to  believe  that  a  nuclear  cloud  has  a 
grossly  nonuniform  structure  with  respect  to  soil  content  and  temperature.  Of 
particular  significance  is  the  sharp  temperature  gradient  known  to  exist  in  the 
cloud.  The  cloud  is  very  hot  in  the  region  of  the  ring  through  the  cloud  about 
which  the  toroidal  circulation  is  centered.  Beyond  this  region  in  all  directions, 
temperature  falls  off  rapidly  with  distance.  One  would  expect  to  find  essentially 
all  of  the  radioactive  device  debris  in  the  hot  region  of  the  cloud.  Also,  within 
this  region  one  would  expect  to  find  whatever  vaporized  soil  material  is  in  the 
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cloud.  As  the  cloud  rises,  this  hot  material  continues  to  expand  into  its  surroundings 
and,  thereby,  to  heat  and  contaminate  additional  soil.  In  light  of  this  description  of 
the  cloud  structure  it  is  apparent  that  actually  it  is  inappropriate  to  use  temperature 
and  soil  concentrations  obtained  by  simple  averaging  over  the  total  cloud  volume. 

Also  of  importance  here  is  the  expectation  that  much  of  the  soil  burden  in  the  hottest 
region  of  the  cloud  will  not  be  in  equilibrium  contact  with  the  hot  gases  there;  in 
other  words,  we  can  expect  a  considerably  lower  average  temperature  to  prevail 
for  soil  than  for  the  gases  in  the  cloud. 

From  Figure  1  it  can  be  seen  that  immediately  subsequent  to  t{  the  cloud 
begins  to  rise  upward  rapidly.  A  strong  updraft  is  produced  in  the  wake  of  the 
cloudk  and  soil  dust  suspended  in  air  is  sucked  upward  in  the  stem  toward  the  cloud 
cap.  The  vortex  circulation  in  the  cloud15’ 16  is  such  as  to  cause  very  high  up¬ 
ward  circulation  velocities  to  occur  around  the  center  axis  of  the  cloud  torus.  This 
tends  to  pull  some  stem  material  into  the  bottom  of  the  cloud  and  hence  to  incor¬ 
porate  this  material  into  the  outer  cloud  layers.  This  means  that  the  cloud  soil 
burden  probably  is  incomplete  at  precisely  tj.  We  can  expect,  however,  that  the 
burden  is  completed  very  shortly  thereafter.  Unfortunately,  there  are  not  avail¬ 
able  quantitative  data  on  internal  cloud  structure  at  times  as  late  as  those  of  inter¬ 
est  here,  nor  do  we  know  how  much,  if  any,  soil  is  incorporated  into  the  cloud 
after  t^. 

In  the  paragraphs  that  follow  and  in  Appendix  B,  we  present  the  results  of 
tentative  estimates  for  the  fraction  of  the  pofl  burden  vaporized  and  the  temperature 
of  the  condensed  phase^soil.  This  work  was  done  after  publication  of  our  initial 
conditions  final  report3  because  additional  information  on  cloud  properties  had  be¬ 
come  available,  as  mentioned  earlier.  Specifically,  we  have  determined  cloud 
altitudes  and  dimensions  as  a  function  of  yield  and  height  of  burst  at  the  initial 
time.  This  allows  energy  balance  calculations  to  be  made  snd  we  hav"  applied 
these  considerations  to  a  determination  of  the  desired  quantities.  We  repeat  that 
the  information  sought  frequently  is  lost  in  the  uncertainties  that  arise  because  of 
the  approximations  used  to  describe  the  energy  partition  model.  The  results  are 
in  a  very  real  sense  arbitrary,  since  by  making  different  basic  assumptions  in 
the  energy  model,  very  different  results  can  ensue.  Furthermore,  we  must  note 
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that  of  all  poaaible  plausible  energy  partition  models,  we  have  had  the  opportunity  to 
investigate  only  one;  fortunately  this  model  does  give  physically  reasonable  vesults. 
Before  proceeding  into  a  discussion  of  the  energy  balance  calculations  let  us  consider 
some  of  the  general  properties  of  the  system. 

Their,  rre  two  categories  of  surface  or  soil  material  to  be  considered:  siliceous 
and  calcareous.  Characteristic  examples  of  siliceous  soil  are  the  common  sandy  and 
clay  soils.  In  the  case  of  sand,  the  chemical  composition  is  almost  pure  silicon 
dioxide  (SK^).  For  the  much  more  abundant  clays,  the  composition  is  more  com¬ 
plex  and  there  are  significant  quantities  of  a  variety  of  metals  such  uu  aluminum 
and  irr*i;  but  even  for  clays,  the  major  chemical  constituent  is  silicon  dioxide.  An 
av^re^e  composition  for  igneous  rock,  from  which  the  common  clay:  'e  derived, 


Constituent 

Percent  Weight 

Constituent 

Percent  Weight 

ao2 

60.18 

CaO 

5.17 

Ai2°3 

15.61 

Na2° 

3.91 

FV>3 

3.14 

*2° 

3.  19 

FeO 

3.88 

TiO„ 

41 

1.06 

MgO 

3.56 

P2°5 

'J.  30 

We  see  that  the  major  constituents  of  clay  soil  (at  high  temperatures)  are  silica  and 

25  2fi 

alumina.  Thermal  properties  for  these  materials  arc:  ’ 


Compound 


ao. 


a*2o3 


Melting  Temperature  Boiling  Temperature 
(°K) _ (°K) _ 

2000  3070 

2323  3773 


Common  examples  of  calcareous  surfer*  materials  are  limestone  and  coral. 
These  substances  are  essentially  caloiu n  carbonate  (CaC03).  However,  at  rela¬ 
tively  low  temperatures,  calcium  carbonate  loses  carbon  dioxide  gas  to  yield  a 
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refractory  residue  of  calcium  oxide  (CaO).  Therefore,  in  considerations  of  cal¬ 
careous  material  at  high  temperatures,  we  use  the  properties  of  calcium  oxide; 
melting  temperature  2853°K  and  boiling  temperature  3123°K. 26 

We  see  then  that  the  minimum  boiling  temperatures  of  siliceous  and  calcareous 
soils  are  approximately  3000°K  and  3100°K,  respectively  Figure  5  shows  that 
cloud  temperature  at  tj  sometimes  is  less  than  the  boiling  temperatures  of  the 

soil  material.  In  such  a  case  we  simply  take  the  fraction  of  soil  in  the  vapor 
phase  to  be  zero. 

In  estimating  limits  on  the  fraction  of  vapor  phase  soil  for  those  cases  where 
the  cloud  temperature  is  above  tha‘  at  the  soil  boiling  temperature,  we  use  the 
following  realizations: 

1.  Temperatures  in  the  cloud  above  the  soil  boiling  temperature 
are  localized  to  th  :  region  of  the  vortex  ring  core. 

2.  This  (thermally)  hot  region  of  the  cloud  contains  virtually  all 
of  the  cloud  radioactivity. 

3.  For  surface  bursty  a  large  portion  of  .alout  is  derived  from 
unvaporized  soil  particles  on  which  radioactive  weapon  debris 
material  condenses. 

4.  The  amount  of  soil  vapor  in  the  cloud  will  be  less  than  the 
mas  i  of  active  fallout  produced  by  the  detonation. 

We  assert  that  item  4  follows  from  items  1,  2,  and  3.  At  the  vaporty.ed  weapon 
debris  material  (which  has  a  negligible  total  mass  compared  with  the  total  soil 
burden  mass)  resides  in  the  hot  core  ring  of  the  cloud  along  with  the  vaporized 
soil.  Both  the  debris  material  and  the  vaporized  soil  subsequently  coudense  on 
solid  or  liquid  soil  particles  to  produce  the  radioactive  portion  of  the  cloud  soil 
burden.  Obviously  the  fallout  material  must  include  a  significant  amount  of 
material  that  was  never  vaporized  in  addition  to  essentially  all  of  the  vaporized 
soil  material.  Therefore,  the  fraction  of  the  total  fallout  mass  that  is  found  to 
be  radioactive  must  exceed  the  fraction  of  vaporized  soil  in  the  c\oud. 
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A  survey  of  the  literature  to  determine  mass  percent  of  active  material 
observed  in  fallout  has  revealed  that  most  of  the  mass  and  most  of  the  activity 
ave  associated  with  particles  in  the  large  particle  size  range  of  the  size-frequency 
distribution.  For  this  reason,  and  since  on  the  average  there  is  only  a  slight 
tendency  for  specific  activity  of  gross  fallout  to  change  with  particle  size,  re 
shall  concentrate  our  attention  here  on  the  larger  particles.  Table  1  i*  a  sum¬ 
mary  of  the  best  data  available  for  large  particles.  We  see  that  an  average  of 
approximately  50%  of  gross  fallout  ma.-s  is  associated  with  active  material.  There¬ 
fore,  on  the  basis  of  the  reasoning  presented  above,  we  can  state  with  a  reasonable 

degree  of  confidence  that,  on  the  average,  less  than  50%  of  the  total  soil  burden  is 
vaporized. 


TABLE  1 


SELECTED  MASS  AND  ACTIVITY  DATA  FOR  LARGE  PARTICLES 


Operation  / 

;==== 

Approx.  Mass 

Approx.  Mass 

/ 

Particle 

of  Total 

of  Total  Fallout 

Source 

/ 

Diameter 

Fallout  in 

Associated  with 

of 

/ 

Range 

Size  Range  ' 

Active  Material 

Data 

y  shot 

W) 

<%) 

j  in  Size  Range 
<%) 

REDWING 

Zuni 

>  420 

75 

75 

Ref.  20 

Fig.  3.53-3.55 
Table  3. 17 

REDWING 

Tewa 

>  210 

90 

75 

Ref.  20 

Fig.  3.57-3.59 
Table  3. 19 

JANGLE 

Uncle 

>  210 

15 

Ref.  22 

Table  IV 
(Station  E) 

SUN  BEAM 
Small  Boy 

200-1000 

58 

27 

Ref.  27  * 

Table  8. 2 

SUNBEAM 
Johnie  Boy 

>  200 

80 

50 

Ref.  28 

Table  3. 17 

AVERAGES 

76 

48 

Only  data  from  collection  stations  in  the  heavy  fallout 
These  were  stations  305,  405,  and  503. 


area  were  considered. 
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To  be  able  to  determine  a  more  definite  value  for  the  fraction  of  the  soil  bur¬ 
den  vaporized,  and  at  the  same  time  estimate  the  temperature  of  the  condensed 
phase  material,  we  have  employed  an  energy  partition  analysis  of  the  detonation. 

The  critical  problem  in  our  analysis  is  to  determine  the  fraction  of  the  total  de¬ 
tonation  energy  release  that  is  consumed  in  heating  and  vaporizing  soil.  Since  we 
already  have  relatively  accurate  knowledge  of  cloud  aif  i*ide  (therefore,  pressure) 
cloud  volume,  temperature  of  tM  cloud  gas,  and  total  soil  burden,  a  knowl  dge  of’ 
energy  available  for  heating  and  vaporizing  soil  allows  estimate,  of  the  temperatures 
of  the  cloud  contents  to  be  made  (see  commentary  on  pp.  8-9  ) 

Estimates  of  energies  available  for  heating  soil*  were  obtained  by  differencing 
the  thermal  energies  radiated  by  air  and  surface  bursts. 29  In  Kef.  29  the  thermal 
energy  radiated  is  given  as  a  function  of  yield  for  both  air  and  surface  bursts.  With 
this  maximum  amount  of  energy  available,  it  is  a  relatively  simple  matter  to  de¬ 
termine  the  desired  quantities  so  that  energy  is  conserve.  The  details  are  given 
in  Appendix  B.  The  results  obtained  are  as  follows: 

1.  1. 5%  of  the  soil  burden  is  vaporized  per  100°  excess  of  the 
initial  temperature  (as  calculated  from  Eqs,  (6),  (7),  and 

(8)  over  3000°K  for  siliceous  soils  and  3100°K  for  calcareous 
soils. 

2.  The  temperature  (°K)  of  the  un vaporized  soil  is  given  by 

T  -  501og1()W  ♦  K00 

where  W  is  the  detonation  energ/  yield  In  kilotons. 


Article  «ze  Frequency  Distribution 


The  prob’mn  o'  determining  pa  tide  size-frequency  distributee  at  t.  is 
essentially  the  problem  of  determining  the  extent  to  which  the  preshot  soil  size- 
frequency  distribution  ha.  b~n  altered  by  vho  proce.se.  of  evaporate  and  growth 
In  these  considerations  we  can  make  a  broad  distinction  between  the  high-yield  and 
low-yield  shots.  We  will  consider  higfa-yield  shots  first.  At  times  prior  to  t  ,  the 
fireball  1.  extremely  hot  (see  Figure  9)  *  „|  because  of  the  stable  nature  of  the  vor¬ 
tex  circulation,  which  during  the  early  *  «e.  of  the  cloud  development  tends  to 


* 

«!] j!  *5!?  eXj"1t  0ther  mean*  for  estimating  the  available 
east  as  plausible  and  can  give  different  resulto. 


ann*  ikni 
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Injure  that  circulating  v  serial  will  remain  contained  within  a  limited  region  of  the 
cloud,  we  can  expect  realdeoce  times  of  entrained  material  In  the  hot  core  to  be 
long  enough  so  that  virtually  all  of  this  material  will  be  completely  vaporized.  Be¬ 
yond  this  region,  condensed  phase  material  should  have  essentially  the  preshot  soil 
size -frequency  distribution.  Subsequent  to  tj,  additional  soil  may  be  entrained  into 
the  cloud;  of  course  this  also  will  have  the  preshot  sell  size -frequency  distribution. 
Subsequent  to  ^  we  can  expect  growth  to  occur  mainly  by  two  mechanisms:  (1)  by 
condensation  of  vapor  onto  toil  particle*,  and  (2)  by  agglomeration  of  melted  or 
partially  melted  soil  particles.  For  high-yield  shots  we  can  expect  both  of  these 
mechanisms  to  yield  significant  growth.  The  condensed  phase  material  present 
as  early  a*  tj,  however,  would  bo  expected  to  have  essentially  the  preshot  soil  size- 
frequency  distribution. 

For  small  detonations  we  note  from  Figure  5  that  the  "average"  temperature 
at  tj  la  below  the  boiling  temperature  of  the  soil  material.  We  must  assume, 
therefore,  that  there  Is  essentially  no  soil  vapor  present  at  tj.  Also  we  note  from 
Figure  3  that  for  the  small  detonations,  tj  is  small;  In  other  words,  the  time  scale 
of  events  is  greatly  compressed  for  small  detonations.  This  means  that  there  is 
much  less  time  available  for  thermal  equilibrium  between  soil  and  cloud  gases  to 
be  reached.  We  hypothesize,  ard  this  la  supported  by  energy  conservation  calcu¬ 
lations,  that  relatively  little  soil  la  ever  vaporized  in  ima)'  detonations;  conse¬ 
quently,  there  cannot  have  been  significant  particle  growth  by  condensation  from 
tbs  vapor.  Furthermore,  since  growth  by  agglomeration  la  a  alow  process,  and  the 
cloud  cooling  rate  la  very  fast  in  relation  to  it  (as  shown  by  cloud  rise  simulation 
calculi  ions),  we  conclude  In  light  of  the  compressed  time  scale  that  growth  by  th’ 
mechanism  also  is  not  significant. 

In  summary,  we  find  that  the  preshot  aoil  size-frequency  distribution  is  an 
adequate  approximation  to  the  aize-f requeue y  distribution  of  condensed  phase 
material  In  the  dead  at  tj.  As  a  corollary  to  this  finding,  we  anticipate  that  sub¬ 
stantial  particle  growth  occurs  only  for  large-yield  detonations. 


29 


COMPUTER  PROGRAM  OUTLINE 


TAILS  I 

a  nummary  or  fnocmam  rvncnom 


PKOCMUJf 

run 

CALIXHO 

mOtMAM 

FUNCTION 

LfNKI 

m»cvm«  uu  data  non  mm  taw.  caul 
oodiwtatioh  mnmtnrrmm,  mrr  m  »  multi 

or  mriiAi.  <x>nmtio»  comfutatiom,  ajb 

TmniATI  MOfdLT  Tt  mW  CALLOW  LOUtJ 

OJ  M KNELT  tTOPFlNO. 

rait 

umu  | 

prroumri  to*  tdii  or  ootul  conormi 

IFtanCATWA. 

mtp 

LfNKI 

tmhwhob  or  fafcm  aod  m 
TAWWinj  ML 

MAM 

LfNKI 

DrTTFMIDP.  toil  1 OAMMO  OF  CLOUD, 

VAK* 

LfNKI 

DrnwMmi  mam  or  »oil  vafo*  m  m  cloud. 

MTM 

LOMU 

5LTJ5?  FA*'ncL*  IHI  WtfmwmON  FANA- 
MVTBM,  MEAN  AND  0TANDAND  DEVIATION. 

IF  KBQinARD. 

Figure  9.  Flow  Chart  of  the  Initial 
Condition*  Module  Computer 
Program  Logic  (LINK1) 


Description 

The  Initial  Conditions  Module 
computer  program  consists  of  a  main 
program  and  five  computation  sub¬ 
routines.  Table  2  summarizes  the  namee> 
and  functions  of  these  programs.  Fig¬ 
ure  9  is  a  general  flow  chart  of  this 
group  of  programs.  Because  of  the 
great  simplicity  of  the  programming, 
more  detailed  flow  charts  have  not  been 
prepared.  Instead,  the  discussion  that 
follows  is  provided  to  supplement  the 
FORTRAN  card  listings.  In  its  Initial 
form,  the  DfcLFIC  system  is  designed 
for  execution  on  the  IBM  7094  computer 
via  the  IBSYS-IBJOB  processor,  and 
the  "overlay"  feature  is  used  to  control 
the  input  sequence  of  major  sections  of 
the  system.  To  facilitate  discussions 
of  the  programs,  the  executive  programs 
of  each  major  section  are  assigned 
names;  the  Initial  Coalitions  Module 
executive  program  is  called  LINK1. 

(See  Figure  3  of  Volume  I  for  a  com¬ 
plete  presentation  of  the  executive  pro¬ 
gram  designations. ) 

Program  Discussion 

The  Initial  Conditions  program  is 
designed  to  serve  two  objectives.  First, 
it  can  be  used  to  merely  compute  and 
print  the  values  of  Initial  conditions 
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parameter*  for  an  unlimited  number  of  input  data  aeta.  Second,  it  can  compute  a 
ain^le  aet  of  initial  parameter  value*  and  then  call  the  next  program  (LINK2)  of  the 
DELFIC  ayatem  *o  that  cicud-rise  computations  can  be  carried  out.  We  shall 
now  refer  to  the  program  lieting  to  aee  how  this  la  accompliahed. 

After  setting  two  tape  name  parameters  and  initializing  the  data  set  counter, 
NN,  the  program  reads  an  arbitrary  72  character  identifier  (card)  to  be  used  for 
Identifying  the  initial -conditions  data  set.  Next,  the  data  set  count  N  is  read  and 
an  overall  title  is  written.  At  statement  number  20,  which  marks  the  beginning  of 
the  processing  loop,  the  particle  size-frequency  distribution  specification  option 
parameter,  IDISTR,  is  read.  This  parameter  determines  whether  a  log-normal 
or  arbitrary  (tabular)  particle  size -frequency  distribution  is  to  be  used.  If  the 
latter  option  is  selected,  the  next  datum  to  be  read  consists  of  the  mid-range  dia¬ 
meter  value  of  the  first  (with  smallest  diameter)  particle  *<ze  range  entry  in  the 
tabular  representation  of  the  particle  size-mass  frequency  distribution.  Next,  a 
data,  set  is  read  that  specifies  yield,  height  (depth)  of  burst,  soil  type,  and  finally 
geometric  mean  and  standard  deviation  for  a  log-normal  distribution,  or  if  a 
tabular  distribution  is  to  be  used,  the  number  of  entries  in  the  table  (the  maximum 
number  of  entries  is  40)  followed  by  the  table  itself  (to  be  described  later).  If  a 
log-normal  distribution  is  to  be  used  but  the  mean  and  standard  deviation  data 
card  fields  are  left  blank,  the  mean  and  standard  deviation  are  supplied  by  sub¬ 
routine  D8TBN.  Then,  the  height  (depth)  of  burst  is  converted  from  meters  to 
feet  and  the  scaled  height  (depth)  of  burst  is  computed  and  stored  in  parameter  Z. 
Next,  tests  are  made  to  see  if  the  model  is  appropriate  for  computing  for  the  speci¬ 
fied  Inputs  If  it  is  not,  a  comment  is  made  at  either  statement  150  or  143  and 
eventually  control  reaches  statemant  171.  At  171,  a  test  of  parameter  N  indicates 
whether  the  run  should  be  continued  in  order  to  compute  initial  conditions  for  other 
input  d'-.ta  or  be  terminated.  The  run  is  terminated  whenever  either  all  data  sets 
have  been  dealt  with  or  when  an  error  or  inapplicability  of  the  model  is  discovered. 

H  thfi  touts  at  statements  63  and  66  indicate  the  model  to  be  applicable,  five 
computatbmal  subroutines  -  TIME,  TElfP,  MASS,  VAPOR,  and  DSTBN  -  are 
called  sequentially.  Note  that  subroutine  DSTBN  is  called  onl>  if  no  preshot 
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particle-size  distribution  was  specified  by  the  user.  *  At  statement  96,  certain 
needed  conversions  are  made  and  then  a  complete  exhaustive  printing  of  inputs  and 
outputs  is  made.  Once  agai"  control  reaches  statement  171  and  if  more  data  sets 
are  to  be  processed,  a  transfer  is  made  back  to  statement  20. 

It  should  be  noted  that  all  input  and  output  parameter  values  e;e  held  in  common 
and,  therefore,  the  communication  between  the  Initial  Conditions  Module  (LINK1) 
and  the  cloud  rise  program  (LINK2)  can  be  carried  out  by  merely  leaving  the  LINKl 
common  intact.  Because  LINKl  executes  so  quickly,  no  storage  of  intermediate 
results  on  tape  is  necessary. 

If  the  user  Intends  to  make  use  of  the  particle  growth  option  of  LINK2,  the 
Cloud  Rise  Module,  he  must  specify  a  particle  size  distribution  for  the  preshot  soil 
over  (or  under)  which  the  detonation  occurs.  If  this  is  not  possible,  it  probably  will 
be  necessary  to  employ  a  particle  size  distribution  for  expected  fallout  n  this  .< 
the  particle  growth  option  should  not  be  exercised. 

As  noted  earlier,  the  user  has  the  option  of  specifying  a  log-normal  r  tor 
the  distribution,  or  he  can  specify  an  arbitrary  distribution  in  tabular  fe.m.t  The 
tabular  distribution  must  be  constructed  and  must  be  input  to  the  program  according 
to  the  following  prescription.  The  tabl«  is  constructed  so  that  each  successive  mid¬ 
range  diameter  and  boundary  diameter  is  times  its  preceding  diameter.  Mid¬ 
range  diameters  are  geometric  means  of  the  boundaries.  In  array  WHY (I)  the  mass 
fraction  for  each  size  range  is  specified  in  sequence  of  increasing  diameter.  There 
can  be  no  more  than  40  entries  in  the  table.  A  parameter  DMIN  also  is  specified 
which  is  defined  as 

DMIN  -  Dj  /  (2  Wl)  , 

where  is  the  upper  boundary  diameter  of  the  first  size  range. 

If  the  user  chooses  to  specify  a  log-normal  distribution,  he  should  specify  the 
geometric  mean  (microns)  and  geometric  standard  deviation  (dimensionless)  of  a 
particle  diameter-number  frequency  distribution  function,  and  optionally,  the  param¬ 
eter  DMIN  (see  the  preceding  paragraph).  In  this  case  DMIN  is  used  only  for 

- 

It  is  strongly  recommended  that  the  user  supply  a  particle  size -frequency  dis¬ 
tribution  for  the  preshot  soil  that  is  appropriate  for  the  geographical  location  of  his 
simulated  detonation. 

^Though  provision  is  Included  for  specification  of  an  arbitrary  particle  size- 
frequency  distribution,  the  user  must  specify  a  log-normal  distribution  if  he  intends 
to  use  the  current  version  of  the  DELFIC  Particle  Activity  Module.  (See  Volume  V 
of  this  documentation.) 
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computing  fallout  from  the  cloud  di*  u;*g  its  rise.  It  Is  not  used  for  setting  up  par¬ 
ticle  inputa  to  the  Transport  Module.  If  DMIN  is  not  specified  in  the  input,  a  value 
la  supplied  by  the  cloud  rise  program  (see  Volume  III  of  this  documentation). 

Operating  Information 

Because  of  the  simplicity  of  this  program,  the  following  tabular  representation 
of  the  input  data  deck  (Table  3)  is  deemed  to  provide  sufficient  operating  instructions. 


TABLE  3 

THE  INPUT  DATA  DECK  FOR  THE  INITIAL  CONDITIONS  MODULE 


Card  Sequence 
and  Number 

Content 

Format 

I 

Arbitrary  72  cnaraoier  identitier 

i4a6 

2 

The  number  of  data  sets  for  which  initial  conditions 
computations  are  desired.  If  the  value  of  this  para¬ 
meter  is  zero  or  negative,  a  single  data  set  will  be 
expected  and  atter  computation  of  its  initial  conditions 
parameters,  a  transfer  will  be  made  to  LINK2. 

15 

3 

Preshct  soil  size -frequency  specification  option  para¬ 
meter,  IDISTR.  A  value  of  0  selects  the  log-normal 
distribution  option,  whereas  a  value  of  1  selects  the 
tabular  distribution  option. 

15 

IF 

EDISTR*0  < 

4 

Yield  in  kilotons,  height  of  burst  in  meters,  soil  type 
(1.0  Indicates  siliceous,  2. 0  indicates  calcareous),  mean 
diameter  of  the  log-normal  preshot  particle  size-number 
frequency  distribution  in  microns,  and  the  standard  devi¬ 
ation  of  this  distribution  (dimensionless).  If  no  values 
are  supplied  for  the  mean  and  standard  deviation,  the 
program  will  supply  them. 

5F10.3 

5 

< 

Mid-range  diameter  of  smallest  particle  size  class  to 
be  used  in  representing  the  log-  normal  distribution  (if 
zero,  a  value  will  be  assigned  by  the  program). 

F10.3 

r4 

Yield  in  kilotons,  height  of  burst  in  meters,  soil  type  indi¬ 
cator  (1.0  indicates  siliceous,  2. 0  indicates  calcareous). 

3F10.3 

IF  j 

IDI8TR-I  < 

5 

Mid-range  diameter  of  the  smallest  particle  size  class  to 
be  used  in  the  tabular  representation  of  the  particle  size- 
mass  frequency  distribution  (see  the  particle  size  distri¬ 
bution  discussion  cm  p.  32). 

F10.  3 

6 

Number  of  entries  in  the  particle  size  frequency  table 
(NDSTR  <  40). 

15 

7,8,.#. 

Particle  size  —  mass  frequency  table  entries  with  five 
entries  per  card.  Each  entry  represents  a  mass  fraction 
as  discussed  on  p.  32. 

5F10.3 

N 

Same  as  card  3  for  the  second  data  set 

This  la  height  above  ground  zero.  For  aubaurface  bursts,  the  depth  of  burst 
is  entered  as  a  negative  number. 
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FORTRAN  LISTINGS 


The  FORTRAN  Hating*  are  included  on  pp.  35  through  44. 

SAMPLE  TEST  PROBLEM  AND  PRINTOUT 

On  pp.  45  through  46  two  test  problem  printout*  are  presented:  (1)  with  the 
log-normal  aoil  particle  size  distribution  option,  aud  (2)  with  a  tabular  soil  particle 
size  distribution.  All  inputs  are  printed  along  with  the  output  and  th»  output  is  self- 
explanatory.  Preparation  at  the  input  cards  is  discussed  5,'  computer  Program 
Outline  section. 


nnnnnnnonn 


6IBFTC  LNKI  LIST .DECK .M96/2 

SUBROUTINE  LINK1 

C  INITIAL  CfNDITItNS  (FIREBALL)  MODULE 
C  TECHNICAL  OPERATIONS  RESEARCH  22  SEPT  1966 
C 


LNKI 

LNKI 

LHK1 

LNKI 

LNKI 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

C 

c 

c 

c 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

L 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

r 


LNK1 

PROGRAM  Ti  DETERNINE  the  INITIAL  CONDITIONS  SPECIFICATIONS  BF  LNKI 
TINE.  TEMPERATURE.  TfTAL  SOIL  NASS.  FRACTION  BF  THE  SBIL  BURDEN  INLNK1 
THE  VAPBR  PHASE.  AND  THE  SIZE  FREQUENCY  DISTRIBUTIBN  BF  THE  LNKI 

CONDENSED  PHASE  SBIL  LNKI 

LNKI 

THE  FIRST  CARD  CONTAINS  ANY  ARBITRARY  ALPHANUMERIC  IDENT IF1CAT10N.LNK1 
THE  SECOND  CARD  BF  THE  DATA  DEC*.  CONTAINS  THE  NUMBER  BF  CASES  TO  LNKI 
BE  RUN.  FORMAT  (15).  LNKI 

THIS  PARAMETER  SHOULD  BE  LEFT  BLANK  IF  THE  USER  WISHES  THE  PR0GRAMLNK1 
TO  CALL  LINK2  AND  SHOULD  BE  GIVEN  SOME  POSITIVE  VALUE  N  IF  LNKI 

THE  USER  WISHES  THE  PROGRAM  TO  STOP  AFTER  COMPUTING  N  SETS  OF  LNKI 
INITIAL  CONDITIONS.  lnki 

LNK1 

other  input  parameters  are  -  test  parameter  iidistr)  to  determine  lnki 

IF  THE  PARTICLE  SIZE  FREQUENCY  DISTRIBUTION  IS  LOG-NORMAL  OR  LNKI 

tabular,  yield  in  kilotons.  height  <depth>  of  burst  in  meters.  lnki 

A  SOIL  TYPE  INDICATOR.  DIAMETER  OF  THE  SMALLEST  PARTICLE  SIZE  LNKI 
(MICRONS).  MEAN  (MICRONS)  AND  STANDARD  DEVIATION  FOR  A  LOG-NORMAL  LNKI 
PARTICLE  SIZE  FREQUENCY  DISTRIBUTION  .  OR  IF  A  LNKI 

TABULAR  DISTRIBUTION  IS  USED.  THE  MEAN  AND  STANDARD  DEVIATION  ARE  LNKI 
DELETED  AND  REPLACED  BY  THE  NUMBER  OF  ENTRIES  IN  THE  TABLE.  IF  A  LNKI 
LOG-NORMAL  DISTRIBUTION  IS  TO  BE  SUPPLIED  BY  THE  PROGRAM.  THE  LNKI 
MEAN  AND  STANDARD  DEVIATION  FIELOS  ARE  LEFT  BLANK.  LNKI 

LNKI 

FOR  UNDERGROUND  BURSTS  INPUT  DEPTH  OF  BURST  AS  A  NEGATIVE  NUMBER  LNKI 

LNKI 

THE  OUTPUT  UNITS  ARE  MASS  IN  KILOGRAMS.  LENGTH  IN  METERS.  TIME  IN  LNKI 
SECONDS »  TEMPERATURE  IN  DEGREES  KELVIN.  YIELD  IN  KILOTONS.  LNKI 

DISTRIBUTION  PARAMETERS  IN  MICRONS  LNKI 

LNKI 

. .  GLOSSARY  . . 

LNKI 

INITIAL  CONDITIONS  IDENTIFICATION  ARRAY  LNKI 

MEAN  DIAMCTER  OF  PARTICLE  SIZE  DISTRIBUTION  (MICRONS)LNKl 
DIAMETER  OF  SMALLEST  PARTICLE  SIZE  riCRONS)  LNKI 

CONTROL  INTEGER  FOR  PARTICLE  SIZE  DISTRIBUTION  LNKI 

0  -  LOGNORMAL  DISTRIBUTION  LNKI 

1  -  TABULAR  DISTRi JUT  I  ON  READ  IN  ON  CARDS  (ARRAY  WHY) LNKI 


DETID(I) 

DIAM 

DMIN 

IDISTR 


IS 


ISIN 

ISOUT 

N 

NDSTR 

NN 

RMIN 

SD 

SSAM 


CONTROL  INTEGER  SfFOFIES  WHETHER  LOGNORMAL  LNKI 
DISTRIBUTION  IS  SPECIFIED  BY  THE  USER  OR  BY  THE  LNKI 
PROGRAM  lnki 
0  -  PROGRAM  SPECIFIED  LOG-NORMAL  DISTRIBUTION  LNKI 
1  -  USE,?  SPECIFIED  LOG-MORMAL  DISTRIBUTION  LNKI 
SYSTEM  INPUT  TAPE  LNKI 
SYSTEM  OUTPUT  TAPE  LNKI 
CONTROL  INTEGER  -  NUMBER  OF  INPUT  BURSTS  LNKI 
LENGTH  OF  ARRAY  WHY  lnki 
INTEGER  -  TESTS  NUMBER  OF  BURSIS  RUN  AGAINST  THE  LNKI 
NUMBER  OF  BURSTS  TO  LL  RUN  LNKI 
RAO I US  OF  SMALLEST  PARTICLE  SIZE  (MICRONS)  LNKI 
STANDARD  DEVIATION  OF  PARTICLE  SIZE  DISTRIBUTION  LNKI 
MASS  OF  CONDENSED  PHASE  MATERIAL  AT  SPECIFICATION  LNKI 
TIME  LNKI 


0 

1 

2 

3 

6 

5 

6 
7 
0 
9 

10 

11 

12 

13 

16 

15 

16 
17 
10 

19 

20 
21 
22 
23 
26 

25 

26 
27 
20 

29 

30 

31 

32 

33 
36 

35 

36 

37 
30 
39 
60 
61 
62 
63 
66 

65 

66 
67 
60 
69 

50 

51 

52 

53 
56 

55 

56 

57 
50 
59 
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VPR 

W 

WHY! I) 


IA2£.*F  INITIAL  CONDITIONS  SPEC '  f .  *  ION 

AVERAGE  TEMPERATURE  OF  GAS  IN  CL»  h 

-JJAGE  TEMPERATURE  if  CONDENSED  Th/'E  MATERIAL  IN 

TEMPORARY  STORAGE 
SOIL  CLASS  INDICATOR 

i*o  for  siliceous 

0.0  FOR  CALCAREOUS 

VAPOR  IN  CLOUD  AT  SPECIFICATION  TIME 
WEAPON  YIELD  <KT> 

c?tST,#M  ,F  T,TAL  articulate  MASS  IN  I-TH 
u^?JirLLSI2!  CLASS*  HAAIMUM  LENGTH  of  ARRAY  -  40 
HEIGHT  OF  BURST  (METERS) 

SCALED  HEIGHT  OF  BURST 


COMMON  /SET1/ 

1  »DETIDI12I »IRISE  •  IEXEC  •  IS1N  •  ISOUT 

?  *  SPAR  .  SSAM  ,  TME  .  TMP1  .  TMP2 

f  T2M  •  U  .  VPR  •  W  »  X  i  2 

5  SPARi0’*  *  SPAA1  *  SPAR2  •  SPAR3 

5  SPAR*  •  RMI.I  ,  SPARS  •  SPAR?  •  SPARS  •  SPAR9 


LNKI  SO 
LNK1  SI 
N  LNKI  62 
LNKI  S3 
LNKI  64 
LNK1  S3 
LNKI  SS 
LNKI  S7 
LNKI  SS 
LNKI  69 
TH  LNKI  70 
LNKi  71 
LNKI  72 
LNKI  73 
LNKI  74 
‘•••LNKi  75 
LNKI  76 
LNKI  77 
•  LNKI  76 
»  LNKI  79 
»  LNKI  SO 
’  LNKI  SI 
LNKI  62 
LNKI  S3 


. . LNK!  84 

FORMAT (12AS I  LNK1  •* 

iSkt'”;?0""*  S,,'C,,IE0  5T"<°"*=  15  WM.lve  HENCE  itSI  “ 

FORMAT  I 5F10. 3 )  LMKl  •• 

FORMAT  ( /  //  23X28H****  INPUT  PARAMETERS  ••••/20X,5HYIELD.4I,X.F1  *  11 

«?«:EE,0HT  *" °'PJH " 

FORMAT  ( 1H+.63X , 9HSILICE0US  )  H**,1 

FORMAT (1H*.63X.10HCALCAREOUS)  ,*?}  ?3 

F3RMATI//20X30HPRE-SH0T  SOIL  PARTICLE  SUE  FREQUENCY  DISTR 'NUTION/LMki  11 

125X32HA  LOG-NORMAL  DISTRIBUTION  WITH  ~/3oUh£aN.31X  ,£l2.^J  11 

2ICR0NS/30X  •  lfHSTANDARD  DEVIATION,  17X, £12.3  /2^  sJh^S^i  li 

3S  DISTRIBUTION  WAS  SPECIFIED  BY)  /23a»34HTHILNK1  97 

FORMAT ( 1H+ ,65X  ,  1 1HTHE  PROGRAM)  ,!!?}  I? 

FORMAT  1 1H+  *S5X ,BHlH£  USER)  TH?, 

FORMAT!  15)  1°<> 

/Smew"?”"'™'  ic*lE°  MP,H  ,F  ls  THE  SCWE  tf  THELNK1  }« 

SCALEO  HEIGHT  OF  DURST  IS  SUCH  THAT  THERE  IS  NOLNki  ioa 
1  *®,L  KASS  tNTRAINED  IN  THE  CLOUD  AND  HENCE  NO  LOCAL  FALLOUT///)  LNKI  lot 
FORMAT 1///23X37H....  INITIAL  CLOUD  PROPERTIES  AT  H  ♦E12«5«14H^SECLNK1  id 
10NDS  ***•/ /20X ,23HAVERAGC  GAS  TEMPERATURE3SX.E12.5.2X.14H0EGRFF«  i  « ! 
2KELVIN//20X.5SHAVERAGE  TEMPERATURE  OF  CONDENSED  PHASE  MATER  I  JfE  IN  tit!  !SI 
3CLBUD,5X,E:2.5.2X.14HDEGR:ES  KELVIN//20X,31)«ASS  OF  VAPORI/Fn  tin. 

4  IN  CLOUD, 30X,E12.3,2X.9HKlL0GRAMS//20X41MiAsS  OF  C*J!SEEJ5C?Hi!J  ^} 

^MATERIAL  IN  CL0LD,2OX,E12.3,2X  ,  9HK I  L0GRAMS//20X,84HP  ARTICLE  SUE  FLNK1  111 
CRiWNcr  DISTRIBUTION  AT  THE  TINE  tf  INITIAL  CMoiVlMS  s«ctn!»5t!2!  “j 

F0RMAT(1H1///S1X,14H*«»»  data  set  12, SH  •*••///) 

FORMAT  I1X.14HLEAVING  LINK  l; 

FORMAT!  1HI///51X19H*  . .  »//12X101HT  HE  n  E  P  A  *  rfEr? 

1  M  e  N  T  OF  0  E  F  E  N  S  E  r*  A  L  LOUT  ?  5  E  D  l  C  ?  I  ? 

?  N  S  Y  S  T  £  M,//31X,19H*  ••••••••  •////43X.36HINI  TIAL  J  f}  2 

3NDITI0NS  (FIREBALL)  MODULE/ //35X, UNPREPARED  BY/oJ Ow{5m&C£S}  III 


36 


i»c»ImEa///*°M,MC  SPEC,ntD  sue  is  .eoative  »ekce  a«i  }f] 

II  F0RMAT  t  /  /2OX5OHPRE-SH0T  SIR  PARTICLE  SUE  FREQUENCY  01  STRIBUt  law/MUr } 
12SX41HA  TABULATED  EMPIRICAL  DISTRIBUTION  WITH  -/30X  1  2U/ .21HPARUCLNK1 
*E  CLASSES/30X.25HMINIMUM  PARTICLE  DIAMETER  »10X»EI2»S>',X  . 7HMILNK 1  \ 
3CR0NS/23X.43HTHIS  DISTRIBUTIIN  WAS  SPECIFIFO  BY  THE  USFRl  . 

19  F.RMAT.25X.37MTHE  DISTRIBUTIIN  I S  L0 G-N0RMAL  SIth -/WX  -LiMfAN.31x2  III 
J * i 5 • 2X * 7HM I CR9NS/ 30X • 1 8H5T ANOARO  DEVI AT10N.17X.E12.*  LNKS  i?g 

191  St!SfiET0UMeTE"*10X*«2.»:«:Si|r“.4N-  ,  LHkI  3? 

ill  ss YuAXiiz  vvvr.r.'i  :?e,same  as  that  ^  ln^j  u? 

\n  ?a5Kr;:sji«;;5itsiis?  sue  fr,^™,,  ts  is; 

C  . . 135 

l  * . * . . . . . l!£|  \£ 

NN-I  LNX1  138 

C  21^  ,N,TUL  C0NOITI0NS  RUN  IDENTIFIER  ,L^|  }” 

READ  USIN.1)IDETID(JJ.J-1,12I  ,L?*1 

C  READ  C0NTR0L  INTEGER  LNK1  141 

READ  USIN#10)N  LHKl  142 


C  WRITE  0VERALL  TITLE 

WRITE  ( IS0UV.16) (DET ID) J) .J-l » 12» 

20  READ! ISIN.10HDISTR 
IFIIDISTRI210.210.211 

C  21C  TEST  F0R  PARTICLE  SUE  FREQUENCY  SPEC1FICATI0N  0PT10N 
710  READ! ISIN.3IW.X.U.DIAM.SD  tuiriCATlfN  0PT10N 

READ(ISIN.3)0MIN 

RMIN«0.5*DMIN 

c  THE  5sERESM,YESAT0I22C  L#G_N0RMAL  DISTRIBUTI0N  SPECIFIED  BY 

IFIDIAMI21.21.22 
7.1  15*0 

60  T0  23 
22  IS-1 

60  T0  23 

211  READ! ISIN.3IW.X.U 

Readiisin.3)dm:n 

RM1N«0.5*DMIN 

READ! ISIN.10INDSTR 

READ! IS  IN .3 ) I WHY ( I ) . I ■! .NDSTR  I 


C  23  C0NVERT  H0B  -  D0B  FR0M  ME 
23  X-X/0.3040 

C  2  IS  THE  SCALED  H0B  -  D0B 
60  2*X/((W)**(l«0/3«4)) 


D0B  FR0M  METERS  T0  FEET 


LNK1  139 
LNK 1  140 
LNK1  141 
LNK1  142 
LNK1  143 
LNK1  144 
LNK1  145 
LNr 1  146 
LNK1  147 
LNK1  148 
LNK1  149 
LNK1  ISO 
LNK1  151 
LNK 1  152 
LNK1  153 
LNK1  154 
LNK1  155 
LNK1  156 
LNK1  157 
LNK1  158 
LNK1  159 
LNK1  160 
LNK1  161 
LNK1  162 
LNK1  163 
LNK1  164 
LNK1  165 
LNK1  166 
LNK1  167 
LNK1  168 
LNK1  169 


T*  SEC  ,F  THE  Mf0El-  15  APPR0PRIATE 

*  r I  a 1 66  »63  »63 

63  I  Ft  2-100.0) 70. 70.150 

66  I F( 2+20.0) 143.70.70 
70  CALL  TIME 
CALL  TEMP 

call  mass 

CALL  VAP0R 
I F  C IDISTR ) 90 .90 .96 


LNK1  170 
LNK1  171 
LNK1  172 
LNK1  173 
LNK1  174 
LNK1  175 
LNK1  176 
LNK 1  1*7 
LNK1  176 
LNK1  179 
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CALL  TEMP 
CALL  MASS 
CALL  VAPOR 
IF! ID1STRI90.90.96 
C 

C  TEST  FOR  ACCEPTABLE  SPEC  IF  ICAT  IONS  OF  PRE-SHOT  PARTICLE  SUE 
C  FREQUENCY  DISTRIBUTIBN. 

90  IFIS0I91 .92  *92 

91  WRITE  ( ISOUT .2 I 
GO  T0  93 

92  2  F  ( D I  AM )  94 , 95 . 96 
94  WRITE  ! IS0UT  *17) 

C  93  SH0ULO  THE  RUN  BE  HALTED.  YES  T0  190 

93  IFINI190.190.170 

C 


95  CALL  DSTBN 

C 

C  CONVERT  H0B  -  D0B  BACK  T0  METERS  FROM  FEET 

94  XaX*0. 3046 

C 

C  CONVERT  VPR  AND  SSAM  FROM  (TRAMS  TO  KILOGRAMS 
VPR-VPR/ 1000.0 

C  DURING  COMPUTATION  SSAM  CONTAINS  THE  VALUE  OF 

C  GAS  ANO  CONDENSED  PHASE  MATERIAL  IN  THE  CLOUD 

S$AM*SSAM/ 1000.0- VPR 

C 

C  WRITE  INITIAL  CONDITIONS  RESULTS 
WRITE  I I SOUT .4 ) W > X 
IFIU-1. 01301 .301.302 

301  WRITE  I  I SOUT .5  I 
GO  TO  305 

302  WRITE  I  I SOUT. 6) 

305  IFI I0ISTR 1310. 310. 311 

310  WRITE! I SOUT .7I0IAM.SD 
IF  IISI102.103.102 

311  WRITE! I SOUT .18 INDSTR.DMIN 
WRITE ! I SOUT .193  I 


THE  T0tAL  MASS  OF 


WRITE!  I SOUT. mi  (  I  .WHY!  I  )  .  I  «1  .NDSTR  ) 
GO  TO  104 

103  WRITE  I  I SOUT .8 } 

GO  TO  106 

102  WRITE  { I SOUT  t 9  I 

106  WRITE! ISOUT  *13) TME.TMP1 .TMP2. VPR, SSAM 
IF! IDISTRI116.116.117 

116  WRITE! ISOUT. 19IDIAM, SO. DM4N 
GO  TO  118 

117  WRITE! ISOUT. 191) 

118  WRITE! ISOUT. 192  I 
GO  TO  171 

143  WRITE  I ISOUT. 11  | 

GO  TO  171 

150  WRITE  ! ISOUT, 12) 


171 

170 

180 


ANATuro  ?rrE2£,?£.??ETH€R  T®  CALL  UN*  2’  RETuRN  TO  COMPUTE 
ANOTHER  SET  OF  INITIAL  CONDITIONS  .  OR  EXIT. 

!F!N-1 1200,200.170 
IFIN-NNU90.190.160 
NN*NNi 1 

WRITE  ! ISOUT. 14INN 
GO  TO  20 


LNK1  180 
LNK1  181 
LNK1  182 
LNK1  183 
LNK1  164 
LNK1  185 
LNK1  18b 
LNK1  167 
LNK1  188 
LNK1  189 
LNK  1  190 
LNK1  191 
LNK1  192 
LNK1  193 
LNK1  194 
LNK1  195 
LNK.1  196 
LNK1  197 
LNK1  198 
LNK1  199 
LNK1  200 
LNK1  201 
LNK1  202 
LNK1  203 
LNK1  204 
LNK1  205 
LNK 1  206 
LNK1  207 
LNK1  208 
LNK1  209 
LNK1  210 
LNK1  211 
LNK1  212 
LNK1  213 
LNK1  214 
LNK 1  215 
LNK1  216 
LNK1  217 
LNK1  218 
LNK1  219 
LNK)  220 
LNK1  221 
LNK1  22..’ 
LNK1  223 
LNK1  224 
INK1  225 
LNK1  226 
LNK1  227 
LNK1  228 
LNK1  229 
LNK1  230 
LNK1  231 
LNK1  232 
LNK 1  233 
LNK1  234 
LNK1  235 
LNK1  236 
LNK1  23T 
LNK1  236 
LNK1  239 
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200  glTjMMUT.l,, 

ItO  STOP 
END 


LNIC1  240 
LNK1  241 
LNK1  242 
INK1  243 
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HBFTC  DST8  LIST 

.DECK  »M9*/2 

DoTB 

0 

sueaeuTiNt  ostbn 

OSTB 

1 

C 

D’TB 

2 

3 

c 

OsTB 

4 

CBMMBN  /SET  It 

D^TU 

i 

1  01  AM 

• 

FIC  . 

IRISE 

•  IEXEC 

• 

1  SIN 

• 

IS2JT 

•  DS  IB 

6 

?  SO 

• 

SPAR  . 

SSAN 

•  THE 

• 

T'lPl 

t 

TMP2 

>  OsTB 

7 

3  T2M 

• 

U  » 

VPR 

•  • 

• 

r 

• 

l 

•  DSTB 

a 

4  »HV 

• 

NDSTR  . 

IDISTR 

•  SPAR  1 

• 

SPAR2 

• 

SPAR  3 

.  OSTB 

9 

5  SPAR* 

t 

RMIN  • 

SPARE 

.  SPAR7 

• 

SPARS 

• 

SPARV 

OsTB 

10 

DIMENSION  FIDI12) 

OS  TO 

11 

OIMENSISN  WMT|*C| 

OSTB 

12 

C 

OSTB 

13 

1' 

15 

c 

OSTB 

16 

DIAM-0.*07 

OSTB 

1  7 

SO-*. 00 

OSTB 

IS 

return 

OsTB 

1» 

END 

OSTB 

20 

I 

!! 


u uu  u u uuv 


SIBFTC  MAS  L 1ST • DECK *M94/2  MAS 

SUBROUTINE  MASS  MAS 

MAS 


MAS 

COMMON  /SET  1 /  MAS 


1 

DIAM 

• 

FID 

9 

IRISE 

• 

IEXEC 

9 

I  SIN 

9 

ISOUT 

9 

MAS 

2 

SD 

9 

SPAR 

• 

SSAM 

• 

TME 

9 

TMP1 

9 

TMP2 

9 

MAS 

9 

T2M 

t 

U 

• 

VPR 

• 

M 

9 

X 

9 

2 

9 

MAS 

4 

WHY 

• 

NOSTR 

• 

IDISTR 

9 

SPAR1 

9 

SPAR2 

9 

SPAR  3 

9 

MAS 

9 

SPAR4 

9 

RMtN 

• 

SPARS 

9 

SPAR7 

9 

SPARS 

9 

SPA79 

MAS 

DIMENSION  FIDI12I  MAS 

DIMENSION  WHVIAOI  MAS 

MAS 


*mas 


. . .  •••••••••••••••••••••• 

MAS 

HOO  OR  DOB  MAS 

IEI K 1 2 3D *240  *240  MAS 

2)0  D-2181.S95  MAS 

0 •-/  has 

R-1.125E402M7.S5E-01 I •O- 1 9 .6E-06 I  * I0*»3.0 1-1 9. 11E-12 1 • IO*»5 . 0 )  MAS 
S0.27E*01*(S.S1E-01 1*0-12. 52E-05 l*<U*«3.0> ♦(! .78E-10) • IQ«*5. 01  MaS 
SSAM-  D»IIW»«»I  J.0/J.AH*(«*«2.0I»S  MAS 

GO  TO  2S0  MAS 

240  £-77. 40685  MAS 

SSAM*  E*IIW)»»IJ.0/J.4» (•11180.0-21**2, 01*060. 0*Z|  MAS 

250  RETURN  MAS 

END  MAS 


0 


10 
11 
Ik 
1  3 

14 

15 

16 
17 
10 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 
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81BFTC  TEN  L1ST.0ECK.M94/2 

SUBRBUT INE  TEMP 


( PM  MAN  /SET1/ 


1 

D'AM 

• 

FID 

• 

IRISE 

f 

1  EXEC 

• 

IS,  4 

9 

1  SBUT 

2 

SO 

• 

SPAR 

9 

SSAM 

• 

TME 

• 

T  MP 1 

• 

TMP2 

3 

T2M 

t 

U 

• 

VPR 

9 

VI 

• 

X 

• 

2 

4 

WHY 

• 

NOSTR 

t 

101 STR 

• 

SPAR! 

• 

SPAR2 

• 

SPAR3 

5 

SPAR4 

• 

RMIN 

• 

SPAR6 

• 

SPAR  7 

• 

SPARS 

9 

SPAR  9 

D I  MENS  I  BN  FIDU2I 
DIMENSION  KHY»40I 


•  •••• 


COMPUTE  VAPOR  TEMPERATURE 
Q»2  *¥•'*!  -  ,03921  I 

A«596'J.»C  C1.145)»*IQ/1B0.M*<  I*  1-0. 03948+0. 02637*0/  1 80.0  I  I 

D**-Cr»447J*(M**l  0.04360)1 
TMP1 »A*I I TME/T2M I **B I +1500.0 

COMPUTE  CONDENSED  PHASE  MATERIAL  TEMPERATURE 

TMP2-50.0*AL0G10IWI+ 1400.0 

RETURN 

END 


TEM 

TEH 

TEM 

T£M 

TEM 

TEM 

TEM 

TEM 

TEM 

TEM 

TEM 

TEM 

TEM 

TEM 

TEM 

TEM 

TEM 

TEM 

TEM 

TEM 

TEM 

TEM 

TEM 

TEM 

TEM 

TEM 

TEM 


0 

1 

2 

9 

4 

5 
4 
7 

A 

9 

10 
11 
12 
l  I 

14 

15 

16 
17 
i  B 

19 

20 
21 
22 

23 

24 

25 

26 
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SIBPTC  tlH  UST.0ECK.M9A/2 
SUtMtUTINE  TIME 

C 

c 


CSMMtN  /SET  1 / 

1  DIAM 

• 

no 

2  »0 

• 

SPAM 

5  T2m 

t 

u 

4  WHY 

t 

NDSTM 

«  SPAMS 

f 

MM  IN 

DIMENSItN  MIDI  12  1 

imise 

S&AM 

VMM 

ID  I  STM 
SPAM* 


1EXEC 

TME 

V 

SPAM1 

SPAMT 


I  SIN 
TMP1 
X 

SPAM? 

SPAMS 


I  SMUT 

TMP2 

l 

SPAM  3 
SPAM? 


DIMENSItN  WHVUO) 


£  •••••••••••••••••••••••••••« •«*•••••••••••**•*••••*••••**• •••*•*•••* 

C 

Q«2*W#* | -*0)921  I 

T?M»0.017*  I  I0.045/0.037I»*<Q/180.I  I  •  I  **•  I  0.  49-<  0.07*0/  1  SO.  I  I  I 

TME-l»4.0«T2MI/n***l0.3l I 

METUMN 

END 


TIM 

0 

TIM 

1 

TIM 

2 

*T|M 

i 

TIM 

4 

TIM 

5 

T  IM 

6 

TIM 

; 

TIM 

8 

TIM 

9 

t:m 

10 

TIM 

11 

TIM 

12 

TIM 

13 

»T  |M 

1  A 

»T  | M 

IS 

TIM 

lb 

TIM 

17 

TIM 

18 

TIM 

19 

TIM 

20 

T  IM 

21 

I 
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C  'MN  LIST .DECK »M9*/2 

VAPR 

0 

SUBROUTINE  VAPRR 

VAPR 

VAPR 

1 

2 

a>> *•«»•••••••••••••••••••••••••«••••••••••• •••*•********• •••••**** ###VaPR 

3 

VAPR 

* 

common  /sen/ 

VAPR 

5 

1  0 1  AM  »  FID 

.  IRISE 

*  ifxec 

•  I  SIN 

t 

I  SRUT 

t 

VAPR 

6 

2  SO  •  SPAR 

•  SSAM 

.  TME 

.  TMP1 

• 

TMP2 

• 

VAPR 

7 

3  T2M  •  U 

•  VPR 

•  w 

•  X 

• 

I 

• 

VAPR 

a 

*  WHY  •  NDSTR 

.  101 STR 

1  SPAR1 

.  SPAR 2 

• 

SPAR3 

• 

VAPR 

*> 

9  SPAR*  •  RMIN 

.  SPARS 

.  SPAR 7 

.  SPARS 

• 

SPAR  9 

VAPR 

10 

DIMENSION  F 101 121 

VAPR 

1 1 

DIME MSI  AN  WHVI40) 

VAPR 

12 

r 

VAPR 

13 

1* 

15 

c 

VAPR 

16 

c 

BRANCH  RN  THE  BASIS  OF 

SRIL  CATEGRRV 

SILICEOUS 

TR 

100. 

VAPR 

17 

c 

CALCAPERUS  TR  200 

VAPR 

IS 

IFIU-1. 01100. 100. 200 

VAPR 

19 

c 

• 

VAPR 

20 

c 

IS  THE  crmputed  vaprr 

TEMPERATURE 

higher 

THAN  THE 

siliceous 

SRIL 

VAPR 

21 

c 

BAILING  TEMPERATURE 

VAPR 

22 

100 

1 F ( TMP1-3UOO  >01120)120 

•  110 

VAPR 

23 

no 

VPR«SSAM*0.00015NITMP1 

-JOOC.OI 

VAPR 

2* 

GR  TR  130 

VAPR 

25 

c 

VAPR 

26 

c 

IS  THE  CRMPUTED  VAPRR 

TEHPFRA  TUBE 

4 
U 1 

3 

X 

Than  the 

calcareous 

SRIL 

VAPR 

27 

c 

BAILING  TEMPERATURE 

VAPR 

29 

200 

IFITMP1-3100. 01120. 120 

.115 

VAPR 

29 

ns 

VPR»SSAM*0.00019«tTMPl 

-3100.0) 

VAPR 

30 

GR  TR  130 

VAPR 

31 

120 

VPR-0.0 

VAPR 

32 

130 

RETUR1 

VAPR 

33 

END 

VAPR 

3* 
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APPENDIX  A 


FORM  OF  THE  TEMPERATURE  EQUATION 

Lot  us  visualize  a  rising  nuclear  cloud  as  a  turbulent  eddy  of  hot  iir  that  is 
be'ng  cooled  by  the  mixing  in  of  entrained  ambient  air.  The  effect  of  mixing  on  the 
average  cloud  temperature  can  be  estimated  from  a  calorimeter  equation:  if  AT  is 
the  temperature  change  due  to  the  entrainment  of  mass  AM, 

MTC  *  AMT  C  -  (M  *AM)(T  t-  AT)  C  . 

p  op  ' '  '  p 

where  M  and  T  are  c’<  ud  mass  and  temperature  prior  to  the  entrainment  of  mass 
AM,  To  is  the  ambient  temperature,  and  is  heat  capacity.  This  gives 

(IT  ,  1  (1M 

dt  ‘  (  o)  M  dt 

The  quantity  (l/M)(dM/dt)  can  be  interpreted  as  being  proportional  to  the  reciprocal 
of  a  mixing  time,  t  ,  so  we  may  write 


dT 

(it 


Jr 


(A  1) 


A  characteristic  scale  length  for  the  cloud  will  be  the  radius  R,  and  a  characteristic 
eddy  velocity  will  be  of  the  order  of  the  cloud-rise  velocity  v.  The  effective  turbu¬ 
lent  diffusion  "constant"  would  then  be  of  the  order  of  D  =  Rv.  From  .veil -known 
solutions  to  the  diffusion  equation,  wo  then  infer  that  a  typical  mixing  time  r  will  be 
of  the  order  R  /I)  -■  R/v,  a  not  unexpected  result.  Equation  (A.  1)  becomes 


dT 

dt 


Q’ 


v 

It 


(A .  2) 


The  essence  of  tills  derivaticn  has  been  presented  in  Ref.  A.  1. 


Our  study  of  cloud-rise  data  obtained  from  movie  films  of  nuclear  test  shotsA'2 
has  led  to  the  relations 

(A.  3)* 
(A.  4) 

where  z  is  cloud  center  altitude. 

Substituting  Eqs.  (A. 3)  and  (A. 4)  into  Eq.  (A.  2),  we  obtain 


and 


R  a 


a  * 
t 


dT 

dt 


n 


T  -  T 
_ o 

t 


(A.  5) 


Upon  integration  of  Eq.  (A.  5),  using  conditions  at  the  second  temperature  maximum 
at  the  lower  limits,  we  have 


T  - 


-  T  \ 

O/ 


1  he  i actor  T2m  -  Tr  is  almost  a  constant  that  we  shall  represent  by  K,  and  we 
arrive  at  the  equation 


In  using  Eq.  (A.  6),  we  should  not  forget  that  it  is  based  on  crude  analogies  and  loose 
dimensional  analyses.  Its  major  virtue  is  that  it  fits  the  form  of  the  observed  data. 
Accordingly,  we  should  not  be  surprised  if  certain  of  the  equation  parameters  turn 
out  to  have  values  different  frcm  those  that  we  would  expect  on  the  basis  of  their 
literal  interpretations  as  implied  in  the  equation  derivation. 


♦ 


Equation  (A.  3)  also  has  been  observed  frjm  experiments 
and  has  been  derived  by  TaylorA  6  &n(j  Batchelor.  A .7 


in  water  tanksA  A  5 
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APPENDIX  B 


CHARACTERIZATION  OF  THE  UNVAPORIZED  PORTION 
OF  THE  CLOUD  SOIL  BURDEN 

In  this  appendix  we  describe  the  results  of  simple  energy  conservation  calcula¬ 
tions  foi  nuclear  detonations  that  are  designed  to  yield  information  on  the  relation 
of  the  temperature  of  unvaporized  soil  at  tj,  the  time  of  initial  conditions  specifica- 
f'on,  to  the  combined  influences  of  detonation  energy  yield,  fraction  oi  this  energy 
remaining  in  the  cloud  at  t.,  height  of  burst,  and  the  fraction  of  the  soil  burden 
vaporized.  An  example  is  presented  of  how  this  information  can  be  used  to  obtain 
estimates  of  the  fraction  of  the  soil  burden  vaporized  and  the  average  temperature 
of  the  unvaporized  portion  of  the  soil  burden. 

Energy  Conservation  Calculations 

Our  energy  conservation  calculations  are  based  in  part  on  the  following 
assumptions: 

1.  Each  kiloton  of  detonation  energy  yield  is  equivalent  to 
4. 2  x  1012  J  (Ref.  B.l,  Table  1.41). 

2.  The  average  temperature  of  the  gaseous  contents  of  the  cloud 
at  t{  is  given  by  Eqs.  (6),  (7),  ard  (8)  of  this  volume. 

3.  The  gaseous  contents  of  the  cloud  obeys  the  ideal  gas  equation 
of  state. 

4.  The  total  mass  of  the  soil  burden  is  given  by  Eq,  (10)  or  Eq. 

(15). 

Also  of  critical  importance  are:  specific  heats  c&  (J/(deg-kg))  for  air  (Pdf. 

B.  2),  where 


946.6  +  0.1971T,  T  <  2300°K 
-3587.5  +  2.125T,  T  >  2300°K 


(B.l) 


55 


and  for  soil  (quartz) 


c 

s 


the  heat  of  vaporization  per  kilogram  of  soil  (SiC>2)  (Ref.  B.4),  Ah  ,  where 

Ahy  =  1.0313  x  107  -  482. 3T  ;  ^  (B-3) 

and  cloud  volume,  V  =  f  .R2H,  where  R  and  H  are  respectively  the  horizontal  and 
vertical  cloud  radii  These  radii  have  been  determined  as  functions  of  time  and 
yield  from  analyses  of  photographic  observations  of  nuclear  test  shot  clouds  (Ref. 

B.  5).  At  tj,  the  cloud  volume  (cubic  meters)  is  given  by 

V  =  1.  78  x  106W0' 928 

’  (B.4) 

where  W  is  the  detonation  yield  in  kilotons.  The  cloud  center  height,  Z,  (for  surface 
and  airbursts  )  (meters)  at  is  given  by 


(Ref.  B.  3),  c  where 
s 


7 

781.  6  -i  0.5712T  -  i^--81.  X-LL- 


T 


T  •  848UK 


(B.2) 


J003.8  +  0.  1351 T, 


T  >  848  K 


Z  =  h  +  108W 


where  h  is  the  height  of  burst  in  meters.  Pressure  at  the  altitude  Z  is  taken  as  tha 
of  the  ARDC  model  atmosphere  (Ref.  B.  6). 


m  t  U  u  the  quartz  data  only.  Specific  heats  for  calcium 

oxide  (CaO)  are  about  20%  less  than  for  Si02  so  that  for  the  high  yield  Pacific  shots 
we  have  overestimated  the  energy  required  to  heat  the  soil  by* this  amount. 

.  n  ^!i0n  ?n!^y  n?1  been  inc!uded  in  the  calculations.  According 

co„CmSed“  tlA,h/ne‘il3^e:°  C°m,“red  ^  ““ 
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If  the  vaporized  fraction  of  the  total  soil  burden  is  specified,  then  the  mass  of 
the  air  in  the  cloud,  ma,  is  given  by 


m  =  M 
a  a 


(B.  5) 


where  Mg  and  Mg  are  the  average  molecular  weights  of  air  and  soil  respectively, 

ms  is  the  ma88  of  soil  gas,  R  is  the  universal  gas  constant,  P  is  pressure,  and  T 
is  the  average  temperature  of  the  cloud  gases.  g 

If,  in  addition,  the  total  energy  available  for  heating  and  vaporizing  the  cloud 
contents  is  known,  then  the  average  temperature  of  unvaporized  soil  is  obtained 
from  the  relation, 


T  =  293  + 
c 


fW  -  Q 


L 


c 

m  c 
s  s 


(B.P) 


where  T,  is  the  average  temperature  of  unvaporized  material,  W  is  the  total  detona¬ 
tion  energy  yield,  f  is  the  fraction  oi  the  total  energy  available  for  heating  toe  cloud 
contents,  mg  is  the  mass  of  unvaporized  soil,  and  Q  is  the  energy  (i.e.  heat 
content)  of  the  cloud  gases  relative  to  the  ambient  condition,  Q  is  given  by 

\  '  (m.c,  +  (Tg  -  293)  *  mf  Ahv  .  m.  7) 


where  Ahv  is  the  heat  of  vaporization  of  the  soil. 

Using  these  relations,  detailed  calculations  were  performed  for  values  of  f 
varying  stepwise  from  0. 1  to  0. 9,  and  for  values  of  m®  representing  fractions  of  the 
total  soil  burden  from  0  to  0. 2,  for  the  rr.nge  of  yields  from  0. 01  KT  to  100  MT  and 
heights  of  burst  from  -20  to  150  scaled  feet  relative  to  mean  sea  level.  The  restric¬ 
tion  that  300  <  Tc  <  3000  was  imposed.  Needless  to  say,  these  calculations  were 
done  on  a  digital  computer. 
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These  calculations  served  to  bracket  the  ranges  of  c»  .ions  within  which 
®nerlty  balances  can  be  attained  Illustrative  graphs  are  “'.own  in  Figure  B.  1  (a-e). 
Values  of  are  given  for  those  points  where  energy  bal.  „es  can  be  obtained. 
Where  energy  cannot  be  conserved,  a  *  Indicates  a  surplus  of  available  energy  and 
a  -  indicates  a  deficit  of  available  energy.  It  Is  apparent  from  these  gtaphs,  that 


0  90  1593  —  2093 —  2623 - + 


+ - + - + - + 

|  A-0 


0.70  II 


l«9 - 1577 — 1041  -1551 - + 


+ - + - + 
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Figure  B.  1.  Illustrative  Graphs  that  Show  Relationship 
Between  Yield  and  Energy  Available  for  Heating  for  Sev¬ 
eral  Heights  of  Burst  (See  text  for  further  explanation.) 
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a  wide  range  of  physically  reasonable  possi . ,<„»  m* 

Accordingly,  It  |»  mandatory  that  a  moans  l,o  found  to  select  appropriate  values  for 
these  parameters  for  specified  yield  and  height  of  burst. 


An  Example  Procedure  for  Specifying  f  and  mB 

- _____ . . .  Q 


As  noted  ptevlously,  the  central  problem  in  the  specification  of  f  and  mB  is  in 

determining  the  fraction  of  the  total  energy  yield  available  for  healing  air  and  soil 

in  the  process  of  easting  about  for  a  reasonable  method  for  determining  this  energy 

ract.on,  we  found  the  following  equations  in  the  19r,r,  edition  „l  the  "Capabilities  of 

Atomic  Weapons"  (Ref.  ».  7)  for  the  thermal  energy  yield  Q..  for  air  and  surface 
bursts: 


^T,a  ~  0-44VV°',)4,  air  bursts 


(B.8) 


^T,  s  0.147W,  surface  bursts 
In  our  example  procedure  we  assume  the  relations 


(B.9) 


fW  =  /  Q  -  Q  \  u 
1  i .  a  ^T,  s  h 


(3.10) 


m  >  0  only  if  T  ^  3000° 


(B.  11)* 


H  is  a  factor  to  account  tor  height  of  burst  that  has  the  form,  of  the  height  of  bu-sl 
sealing  factors  used  for  the  soli  burden  calculations  (Eqs.  (10|  and  (15)): 


H  =  K'WDtz) 
(1X2.  7)2f32.  7) 


A  <  0 


„  .  m  -  A)  (360  ♦  At  t  /[  >  o 
(180)"(360) 


(B. 12) 


>  0  orJyTT8>f°310oeCeOU8  60'‘S'  F°r  Calcare0,JS  »»««• 
8  ff 


we  assume 
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v  here  A  la  given  by  the  height  of  buret  (ft)  divided  by  Wl/:|  4  a„(l  z  _A  F(>r 
Tg  <  3000°,^  we  found  that  energy  balances  are  obtained  for  all  pertinent  yields. 
Above  .1000°,  we  computed  values  of  T,  for  values  o.  m«  obtained  from  the  relation 


where  ms  is  the  total  soli  burden  and  p,  the  fraction  of  the  soil  burden  vaporised 
per  100  excess  of  the  Initial  temperature  over  3000°,  ranged  from  0  to  0,017. 

The  solid  curve  of  Figure  B.  2  shows  unvaporized  soil  temperatures  computed 
as  Indicated  above.  This  curve  remains  substantially  the  same  for  any  height  of 
burst.  For  yields  greater  than  about  103  KT.  vaporized  soil  is  predicted  to  remain 
in  the  cloud  :i  t|.  Beyond  this  yield  we  have  plotted  temperature  curves  for 

several  vain,  p  from  Eq.  (B.  13).  According  to  these  curves,  temperature  of 

unvaporlzed  ncreases  as  p  Increases.  This  Is  an  artifact  of  the  particular 

mode  of  setting  up  the  energy  partitioning  and  Is  illustrative  of  the  sensitivity  of  the 
calculation  to  even  minor  variations  In  such  procedures.  Since  the  overall  sp-ead 

of  temperature  Is  only  300°  and.  In  addition,  the  model  design  Is  rather  arbitrary, 
we  need  not  be  particularly  concerned  by  it. 


The  dashed  curve  of  Figure  B.  2  shows  the  temperature  z-elationshtp  chosen 
(rather  arbitrarily)  to  provide  values  of  T.  for  fallout  prediction  calculations;  it  is 
described  by  the  relation 


Tc  -  50  log1()  W  +1400  .  (B.  14) 

This  expression  yields  satisfactory  energy  balances  except  tor  low  yield  subsurface 
bursts.  For  example,  at  W  =  10"2KT  and  A=  -20,  f  =  1. 16.  That  the  method 
should  fall  for  these  conditions  Is  not  at  all  surprising  since  our  gas  temperature 
values  for  surface  bursts  are  quite  suspect  at  the  extreme  low  end  of  the  yield  range 
and  they  scarcely  apply  at  all  to  subsurface  bursts. 
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The  value  of  u.015  for  the  constant  p  to  be  used  in  fallout  prediction  calcula¬ 
tions  (see  item  1  on  p.  28)  was  chosen  by  the  following  considerations: 

1  •  For  P  as  great  as  0. 017  at  W  =  105  KT,  the  vaporized  soil 
mass  already  represents  as  high  as  28%  of  the  total  soil 
burden  (at  A  ■  150).  Thus,  there  is  little  incentive  to  con¬ 
sider  a  larger  p. 
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2.  According  to  Figure  B.  2,  the  temperature  curve  begins  to 
level  off  at  yield  values  less  than  those  for  which  vaporization 
occurs.  A  value  for  p  of  around  0. 015  seems  to  provide  a 
reasonable  extrapolation  of  this  tendency. 

3.  0. 015  is  a  convenient  number. 

REFERENCES 

B.  I  S.  Glasstone,  editor,  "The  Effect,  of  Nuclear  Weapon.,  ■  USAEC 
(1962,  revised  1964). 

B.  2  j  HUeenrath,  et  ,1. ,  -Table,  of  Hemal  Properties  of  Gases,  • 
NBS  Circular  564  (1  November  1955),  Table  2.3. 

B.3  K.K.  Kelly,  "Contributions  to  the  Data  on  Theoretical  Metallurgy, 
Xm.  High-Temperature  Heat-Content,  Heat-Capacity,  and 

Entropy  Data  for  the  Elements  and  Inorganic  Compounds, »  Bureau 
of  Mines  Bulletin  584  (1960). 

B.4  H.  L.  Schick,  "A  Thermodynamic  Analysis  of  the  High -Temperature 
Vaporization  Properties  of  Silica,  •  Chemical  Reviews  60,  331 
(I960).  ~ 

B.5  H.G.  Nonnsut  and  S.  Woolf,  -stadia,  of  Nuclear  cloud  Rl.e  and 

Growth,  •  Technical  C^eratlon.  Research,  TO-B  68-9  (2  February 
1966).  Secret-F. R.  D. 

B  *  8  L  editor,  -Handbook  of  Geophyalca  and  Space  Enrtren- 

(Air  rente  Cambrtdje  Research  Laboratortoa.  19«5; 
McGraw-Hill  Book  Company,  1966). 

B.  7  "Capabilities  of  Atomic  Weapons,  •  Dept,  of  Army  Tech.  Manual, 

TM  28-100  (1955).  Secret-R.D. 

B.8  j.d.  Coughlin.  "Contributions  to  the  Data  on  Theoretical 

M*rtaUur*y-  3®*  «nd  Free  Energies  of  Formation  of  In¬ 

organic  Oxides,  •  Bureau  of  Mines  Bulletin  542  (1954). 


64 


WIIIIIIM 


Security  Claasilieariea 

OOCUMCNT  CONTROL  DATA  •  R4D 

*•  OWOINAIINO  ACTIVITY  /Corporate  a*jAo  J  _ _ _ _ _ _ _ _ _ _ 

Technical  Operations  Research,  Unclassified LASSFCAYION 

Burlington,  Mass.  -yr-oSouP - - - 

~i  nnpont  titlh - - — - - — - - - — _ _ _ 

Department  of  Defense  Land  Fallout  Prediction  System 
Volume  H  -  Initial  Conditions 


PONT  fCC'JKITV  C 

Unclassified 


*•  OMCNwrivt  NOTH  (Typ,  ofnport  md  Muaiva  ,’a ttt) - “ 

— Final  flsport _ _ _ 

AUTHOP^S^  fLttti  HMfv  /i/lt  XMMf ,,  ^SSSJ* 

H.  G.  Norment,  W.  Y.  G.  Ing,  and  J.  Zuckerman 


T0T*t  NO.  OP  PACKS  ~tK  NO.  OP  ACPI 

- JLZ _  44 


«•"  KKPONT  OAT |  "  - - 

- 30  September  19re _ 

«A  COMTNACt  ON  MANT  NO. 

DAl8-035-AJfC-737(A) 

L  PHOJtCT  I  JO. 

C.  TAM 

DASA  Subtask  A7a/10. 058 

d. 


».  AVAIL  ANILITY/ LIMITATION  NOT  I  CM - - - 1 - - - _ 


M.  SUPPltMKhTANV  NOT  l* 


pnsopino  BiCTTZRT  XZTivTl 


Defense  Atomic  Support  Agency 


t  JAM  M 


UnclsNilfled 


Unclassified 


Security  Classification 


14. 


key  uoros 


Initial  conditions 
Fallout 

Nuclear  weapons  effects 
DELF*C 


LINK  A 


ROLE 


INSTRUCTIONS 


WT 


LINK  0 


ROLE 


WT 


LINK  C 


ROLE 


WT 


of  theRcol^T.,NC  A£nV,TY;  ,he  •"<*  address 

nLuLl  '  r“bcLolUr*c*°f’  r«n‘ee.  Department  of 

the  ^  °rg‘n,“,i0“  <'°n>°ra»  culMor) 

CLASSIFICATION:  Enter  the  ove< 
?.«  1  c“ri*r  -laaaification  of  the  report,  indicate  whether 
Restricted  Data"  is  inched.  Making  i.  to  b,  in  accord- 
aace  with  appropriate  security  regulations. 

DtecURy?U5^AM^iCJTTV,*din^  *!  ,Pecified  DcO 
EW‘tk!  10  cl. Ara?ld  Forces  Industrial  Manual. 
l"®  r«»P  Bomber.  Alao.  when  applicable,  show  that 

m  itthirT^!"*"  h*Ve  b"*"  U*,d  for  ^OUf>  3  ■»«1  Croup  4 
capltal^u JrsT^nai. lh*  c®“Pt,«‘«  «P«t  title  in  all 

irtiLsssi 

t^fisi?wi5S£;‘;vsr.ijtL.**  ■"  c,r“‘u  * 

r»»<2tE^:?IP^iy,LN0TES  11  ■PProPri«**.  enter  the  type  of 
:., *"*?***!  eea^ary,  annual,  or  final 

cowed!  1  wfc*n  ■  •peeitic  reporting  period  is 

fis  “““H01  •*hor<o)  aa  shown  on 
Enter  last  name,  first  name,  middle  initial. 


10.  AVAILABILITY  LIMITATION  NOTICFS-  rn„, 


(1) 

(2) 

(3) 


S3?.""’  «<  ">'• 

‘L\  S.  Government  agencies  may  obtain  cooiea  of 

4k"  ddc 


s.  AUTHOR(S): 

or  ^tha  report 


T*^ch  °*  ee^lee-  The  name  of 
the  principal  author  is  aa  absolute  minimum  requirement. 


«.  REPORT  DATE: 
■oath.  year, 
ea  the  report 


^e1*  of  the  report  as  day, 

r->  ryt’U  amra  th - ^  7 

data  ef  publication. 


7*  NUMBER  OF  REFERENCES, 

referenc  e  cited  ia  the  report. 


Eater  the  total 

RffW  NVMBERtSl:  Eatm  the  otfi- 

^rr,u!i£'  ttIL^uIStj1^ 

•rtf  A,  gpeeserA  a&o  eater  this  aamherfsl.  * 


u’ 

<5)  fi>Jlnnrribu8i0n  k0,Mhl*  rePor‘  «  controlled.  Ouali-’ 
tied  DDC  users  shall  request  through 

H,  S"U»£HT*.Y  NOTES:  l„  („ 

12.  SPONSORING  MILITARY  ACTIVITY-  Enter  th*  » 

**  /.’.rarr.'js'a  fe^^ajrrassV- 

porta  be  2SLlfZld"bF  “k ' '  ,h*  *bW  of  d—ified  r^ 


How. 


-  «r  —  iwv  iv  .AD  WOniSi 

defected  so  that  ao *e7J *wd«  «■*.  be 


-selected  so  .■ 
flora,  such  as 


that  no  aecSriS  -i  re»r  *  *e7  ■■*«  be 

_ ^-Tr"y  flneetltcatwn  is  required.  Ideal  i- 

’my  n.oiect  cod?  g**t*"*ton,  trade  .name,  mili- 

,-r  SMi 


aired.  Ideati- 


Unclassified 

Security  Classification 


